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SOME SUGGESTIONS FOR A NEW THEORY 
OF CELL DIVISION 


N. RASHEVSKY. 


COMMITTEE ON MATHEMATICAL BIOLOGY 
THE UNIVERSITY OF CHICAGO 


A number of shortcomings of the diffusion drag force theory of cell division are 
pointed out, especially in regard to the difficulties of interpreting mitotic phenomena. 
As an alternative a different theory is outlined, based on the emphasis of colloidal 
phenomena in the nucleus. The nucleus is assumed to contain very long molecules, 
which are disoriented in the resting stage but which become oriented prior and during 
the mitosis, forming the mitotic spindle. Such an oriented structure is of the nature of 
a tactoid. The molecular mechanism of the formation of the spindle is outlined, and 
some equations are derived which suggest possible quantitative verification. The 
interplay of the molecular forces is assumed also to be responsible for the splitting of 
the chromosomes. Arrangement of chromosomes in the equatorial plate and their 
movement to the poles are discussed, and equations describing those phenomena are 
outlined. 


The only quantitative theory of cell division which has been developed 
somewhat systematically is the diffusion drag force theory (Rashevsky, 
1948). That theory has some important achievements to its credit. The 
theoretical estimations of average cell sizes agree with actual values (loc. 
cit., chap. xi). The predicted correlation between cell size and rate of 
metabolism (Joc. cit., chap. xi) also leads to a fair quantitative agreement. 
Light is shed on the experimentally observed correlation between gly- 
colitic coefficient and rates of tissue growth (Joc. cit., chap. xv). The most 
striking success of the theory, perhaps, is exhibited by the remarkable 
agreement between the theoretically derived elongation curves of freely 
dividing spherical cells (Joc. cit., p. 164; also Landahl, 1942b) and the 
curves actually observed by R. Buchsbaum and R. R. Williamson (1943). 
A similar agreement is found for the rate of constriction of the “neck” of 
a freely dividing cell, derived theoretically by G. Young WaG. Ctl Ds 148). 

Against these advantages, a number of shortcomings must be listed. 
The principal one is probably the failure of the theory to account for 
phenomena of mitosis. Another shortcoming is the neglect of the role 
played by colloidal phenomena during cell division. mae 

In the early development of the theory it was argued that since mitosis, 
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important as it is, is not a universal accompaniment of cell division, there- 
fore a theory of cell division need not primarily be concerned with a theory 
of mitosis. A theory of mitosis could be developed as a further elaboration 
of the fundamental assumptions, by addition of new assumptions, if neces- 
sary. Such an elaboration of the basic theory was attempted by 
N. Rashevsky (Joc. cit., chap. xvi), without much success. The assump- 
tions made in loc. cit. about the attachments of the fibers to the particles 
seem rather artificial. The conclusion that H. N. Barber’s work (1939) 
supports the theory that chromosome movements are produced by con- 
traction of quasi-elastic fibers (Rashevsky,1941) was shown by F.Schrader 
(1944) to be wrong. Moreover, more and more evidence is accumulating 
showing that there are no continuous fibers in the mitotic spindle. Es- 
pecially convincing seems to be the work of G. Ostergren (1949) on Luzula. 
His conclusions are that the spindle is of a nature of a tactoid or liquid 
crystal, a suggestion also made by J. D. Bernal (1940) and other authors. 

It must be remarked that all possibilities accounting for phenomena of 
mitosis on the basis of the diffusion drag force theory have not been ex- 
hausted. The complex interplay of such forces may conceivably produce, 
without the presence of any fibers, the characteristic stages of mitosis. A 
theory developed along such lines would, however, seem to require a large 
number of artificial assumptions ad hoc. 

The heaviest blow, however, to the line of thought which led to the dif- 
fusion drag force theory of cell division is dealt by the recent work of 
Edward D. DeLamater and his associates (DeLamater, 1952; DeLamater 
and Hunter, 1951; DeLamater and Mudd, 1951, Mudd, Winterscheid, 
DeLamater, and Henderson, 1951), who report observations on mitosis in 
bacteria. If, as these observations may indicate, mitosis is a feature of 
every cell division, then an adequate theory of cell division should account 
for mitosis in a natural way, and not through additional elaboration. 

The other shortcoming mentioned above, namely, the neglect of col- 
loidal changes, is also serious. A number of workers, in particular D. A. 
Marsland and his coworkers (Brown, 1934; Marsland and Brown, 1936, 
1942; Marsland, 1938, 1939; Pease and Marsland, 1939; Pease, 1941; 
Marsland, 1942; Rugh and Marsland, 1943; Marsland, 1948; Marsland 
and Jaffee, 1949), have established a close relation between cell division 
and the gelation of cell colloids. Such a general relation may be rather 
naturally explained in terms of the diffusion drag force theory. The ex- 
pression for the dividing forces [loc. cit., p. 143, eq. (12)] contains the dif- 
fusion coefficient D in the denominator. A gelation results in a decrease of 
D and, therefore, in an increase of the dividing forces. However the details 
of various observations, notably those of D. C. Pease (1941), on the re- 
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appearance of asters after reduction of pressure are hard to reconcile with 
the simple diffusion drag force theory. 

In some cases effects of drugs also present difficulties. While the inhibi- 
tory effect of some drugs on cell division may be explained by their in- 
hibitory effects on metabolism, in some cases the effects strongly indicate 
the involvement of colloidal mechanisms. We may mention here the ob- 
servation of B. Osogoe (1941) that potassium ions result not only in an 
inhibition of division but in an almost complete reversal of a nearly com- 
pleted process of division. ae surely would not totally inhibit the 
metabolic processes. 

Methodologically the diffusion drag force theory-has the following ad- 
vantage. The diffusion drag forces are present in any metabolizing system, 
hence in any cell which manifests life. Since, as the theory shows, those 
forces are sufficient to produce cell division, the theory does represent an 
economy of thought because it does not introduce any extra assumptions. 
In view of what was said above, it seems that the diffusion drag force 
theory of cell division has outlived its usefulness. This conclusion should 
in no way deprecate the possible importance of diffusion drag forces in 
biology (cf. Hearon, 1949a, b; 1950a, b, c). 

Of other theories, the older work of Rashevsky (loc. cit., chap. xiii), 
which was based on the concept of variable surface tension, may be men- 
tioned. This concept did seem to have some experimental justification 
(Spek, 1918). However, though a general theory of division can be devel- 
oped on this basis and some features of elongation during cell division can 
be described mathematically (Shimbel, 1946), the theory always contains 
some indefinite parameters, and no prediction as to size of cells, etc., can 
be made. 

More recently Rashevsky suggested a modification of the diffusion drag 
force theory by considering simple osmotic pressures, produced by dif- 
fusion (Rashevsky, 1949a). This approach has the advantage of leading to 
a theory of segmentation, a phenomenon in which no marked elongation 
is involved and which offers some difficulties to the diffusion drag force 
theory. However the objections which hold against the diffusion drag 
force theory also hold for this theory. 

In another paper Rashevsky (1948b) investigated the possibility that 
elastic stresses, produced by gelation, could result in the division of a cell. 
From considerations of the theory of plasticity, through application of 
Betti’s theorem, he arrives at the conclusion that such elastic forces could’ 
not result in a continuous elongation and eventual division of the system. 
Recently, however, Irwin Isenberg (1952) has shown that Rashevsky’s 
conclusion holds only for the case of a homogeneous system with constant 
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viscosity, for which the classical Betti’s theorem has been derived. By 
generalizing Betti’s theorem to heterogeneous systems, Isenberg shows 
that elongation and division due to elastic stresses may be possible. This 
important work of Isenberg paves the way to a mathematical theory of 
cell division based on elastic stresses in colloidal systems. 

A large amount of qualitative theorizing has been done concerning elec- 
tric charges as possible sources of cell division. The psychological back- 
ground for this theorizing may, perhaps, be found in the rather superficial 
analogy between some mitotic figures and the models of lines of forces 
between electrical charges. The main difficulty with the theory of electric 
forces is that the cell as a whole is electrically neutral. Because of the 
electrical double-layer nature of the cellular charges, rather strong fields 
may exist locally within the cell. But the gross effect on the cell as a whole 
cannot be expected to be sufficiently strong. An estimation of electric 
forces due to metabolic processes which involve diffusion of ions shows 
that they are much weaker than the diffusion drag forces (Joc. cit., chap. 
xili). 

In the following we shall very briefly outline a possible new approach 
to the theory of cell division, which emphasizes both the phenomena of 
mitosis and the role of colloidal structures. 

It has been long known that a cell, deprived of cytoplasm but with 
nucleus intact, will regenerate the cytoplasm; but an ennucleated cell does 
not regenerate the nucleus. This suggests that the nucleus is the basic 
region for biological synthesis, the cytoplasm being formed from the 
nucleus as a sort of by-product or, rather, second step in this process of 
synthesis. The work of T. Caspersson (1950) gives rather direct support 
to this concept. We shall therefore consider the nucleus as the basic eae 
ment of the cell. 

In line with existing evidence (Caspersson, Joc. cit.) we shall consider 
that the nucleus contains long chain molecules P of proteins or other sub- 
stances. We assume that due to intermolecular forces, which may be either 
of the nature of electric dipoles or of ionic cloud origin (Levine, 1939a, b), 
these molecular chains P all tend to orient themselves in one direction. 
Complete orientation corresponds to a minimum of the potential energy 
én. This orientation is opposed by the thermal energy. If m is the number 
of degrees of freedom of a molecule P, then, in order that the orientation 

_be sufficiently pronounced, we must have 


ni 
| em | 2>—— (1) 


where & is the Boltzmann constant and T the absolute temperature. 
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The minimum energy én, whatever its precise nature, is in general a 
function ¢,(¢) of the concentration c of different substances present in 
the nucleus. Ions, especially hydrogen ions, are particularly likely to 
affect ém. If acid substances are produced in the nucleus and if ¢ is the 
concentration of the hydrogen ions, then, for a constant rate c of produc- 
tion, we have 


2 
c= og +28, (2) 


where 7 is the average linear dimension of the nucleus, D the diffusion 
coefficient, and co the concentration outside of the nucleus. 
If |«n| increases with c, so that 


and if c* denotes the concentration above which inequality (1) is valid, 
then there exists a critical size r* above which the P molecules will be 
oriented. That critical size is determined by 


ae oe 
Cot D == Gs (4) 
From (4) we find 
(Gee aT) 
Ne oe (5) 


If we consider as plausible values c* — cy ~ 10-* gm. cm.-*, g ~ 10° 
gm. cm.— sec.—!, and D ~ 10~7 cm? sec.—1, we find r* ~ 10-* cm. Thus 
the orientation will occur in regions which are of the order of magnitude of 
actual nuclei. In developing the theory it will be necessary to make specific 
physical assumptions about the nature of the intermolecular forces, and 
to compute e,, as a function of, let us say, the concentration c of some 
biologically important ions. Then, from the knowledge of the rates of re- 
actions by which those ions are produced and of the diffusion coefficient, a 
more refined equation than (5) should be derived. 

In an anisotropic structure, consisting of oriented long-chain molecules 
(Fig. 1), the surface tension at the “poles” p is in general different from 
that at the “equator” e. We shall assume that the equatorial surface ten- 
sion y. is greater than the polar yp. This assumption implies an inequality 
between the intermolecular forces in the longitudinal and transverse di- 
rections. Our assumption means that it takes less work to remove a mole- 
cule from the equator to infinity than from a pole to infinity. If f(r) is 
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the force acting in the first case and f,(7) the one acting in the second, and 
if d, and d, are the distances of closest approach in the two cases, then we 


must have 


‘ ) Ci aPS ee (6 
i fp(rydr> fi f(r) dr ) 


Pp é@ 


Again, the expressions f.(7) and f,(r) will be determined by the nature of 
the forces. Inequality (6) actually implies some inequality or inequalities 
between the parameters which characterize the mechanism of the inter- 
molecular forces. It will be the task of the theory to investigate under what 
conditions (6) holds, if plausible assumptions are made about the inter- 
molecular forces (e.g., dipoles or ionic clouds). 


FIcureE 1 


We shall consider the chromosomes as bundles of parallel molecular 
chains which multiply and thus increase in thickness. On the side of each 
chromosome there is a molecular group C which has a strong affinity for 
the end of a P molecule. This affinity may be due to a polar character of 
both the group and of P, or it may be due to other factors. The force of 
attraction between C and P is assumed to be greater than the force of 
attraction which holds the individual chains of the chromosome together. 
If ¢, is the minimum value of the binding energy between the chains which 
constitute the chromosome and ¢, the minimum binding energy between 
C and P, then we assume that 


65 Gn (ep) 


The assumption expressed by inequality (6), namely, y. > Y», implies 
that the chromosomes will be arrayed in the equatorial plane along the 
surface. 

If the oriented nuclear structure has a gel-like consistency, then internal 
stresses will develop in it and it may either be deformed, or simply break 
in the region of greatest stress. If the gelation and the orientation are re- 
lated, then the region of maximum stress will be the equatorial region, 
because the concentration c, which is responsible for the orientation and 
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gelation, is largest there. If the stress T(c) at the equator is greater than 
the longitudinal tensile strength S of P, 


Pgye Se (8) 


then the structure simply cracks at the equator. That breaking will partly 
involve the separation of P molecules and partly the separation of struc- 
tures like those illustrated in Figure 2, where Ch denotes schematically the 
component chains of a chromosome viewed “end on.” Because of inequal- 
ity (7), the break will occur in this case not between C and P, but between 
the chains Ch. Hence the chromosomes will be split longitudinally, each 
half belonging to the corresponding “‘half-spindle.” 


FIGURE 2 


The breaking of the spindle results in an increase of its specific surface, 
and therefore in a decrease of c, which may drop below c*. When this 
happens, a disorientation occurs in each half-spindle. This disorientation 
also results in a gradual liquefaction of the structure, which is accom- 
panied by a reduction of volume. (A reduction of volume during solation is 
the rule for cell colloids as evidenced, for example, by the experiments 


SS e e== 
= So 


Ficure 3 


mentioned above of Marsland.) This reduction in volume of the two halves 
will in itself cause a slight initial movement apart of the chromosomes. 

Further movement of the chromosomes may be viewed as due to the 
fact that they are now arrayed at the “poles” of the half-spindles (Fig. 3). 
Since the structures are not yet completely disoriented we still shall have, 
for each half spindle, y. > yp, though, as the disorientation proceeds, 
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Ye — Yp tends to zero. Very roughly, the force pulling each chromosome 
toward the equator of the half-spindle is given by 


(ee (9) 


where / is the length of the half-spindle, and s the area of the cross- 
section of the chromosome. 

The quantity y.— vp is a function of disorientation, and therefore a 
function of time. Its theoretical determination should proceed from con- 
siderations of the “relaxation time’”’ of the oriented structure. This can 
again be roughly estimated by considering the angular velocity of thermal 
rotation of a P molecule. Thus the expression for f in equation (9) will 
involve the size and mass of the P molecules. Once f is known as a function 
of time, the speed of the moving chromosomes can be determined approxi- 


igh Caves —~ Se 
7 —j a 
41} \\ —— ae 
I yA\\ 
a b 
FIGURE 4 


mately from considerations of viscous resistance. The viscosity will in this 
case also decrease as the orientation decreases. The general picture 
adopted here is not unlike that suggested by Bernal (1940), except that 
we do not have to make the additional assumption which Bernal makes, 
namely, that the chromosomes produce a special liquefying substance. 
Moreover, our formulation is somewhat more quantitative than Bernal’s. 

The disappearance and reappearance of the nuclear membrane may be 
connected with the orientation of the P molecules. When the P molecules 
are not oriented, we may have a structure like the one shown in Figure 4a, 
where the shorter molecules, forming the membrane, may be oriented on 
the surface of an isotropic liquid. The radial orientation of the membrane 
molecules may become impossible when the P molecules are oriented (Fig. 
4b), because the whole structure of the field of force changes. 

The important feature of the above suggestions is that they assume five 
imequalities (1), (3), (6), (7), and (8). If we know the range of variation of 
any two quantities « and y, know the mean values and 4, and the stand- 
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ard deviations o, and o, of their distributions, we can evaluate the prob- 
ability of the inequality 
Lee (10) 


If pi(x) and p.(y) are the distribution functions of x and y, then the 
probability of (10) is given by 


plea) =f r(x} fpa(y) dy} ae. (11) 


The range of variations and the distribution of the quantities involved 
in our inequalities (1), (3), (6), (7), and (8) can in principle be determined 
empirically. Therefore, in principle at least, we can calculate the probabil- 
ity of a cell forming by chance in such a way that it will multiply in- 
definitely by division. 

If inequality (1) does not hold, the cell never divides. It will reach 
asymptotically a size limited by the penetration of metabolites. 

If inequality (3) does not hold, then a cell above a certain size will not 
divide. Below a certain size it may divide into the smallest fragments. 

If inequality (6) does not hold, there will be no mitosis. If all chromo- 
somes are essential to the life of the cell, a cell dividing with unequal dis- 
tribution of chromosomes will not survive. The probability of such patho- 
logical divisions may be calculated and compared with observed inci- 
dences, if any. 

Liquefaction or disorientation of the P molecules at early metaphase, 
before the chromosomes are pulled apart, will arrest mitosis. Presence of 
small bodies with strong affinity for the ends of P molecules may result in 
aster-like structures. The mechanism implied in inequality (7) will result 
in the attraction of chromosomes by re-formed “asters,” a phenomenon 
observed by Pease (1941). 

If inequality (1) does not hold, but chromosomes grow without the con- 
stituent chains being separated, giant chromosomes would be formed. It 
would be interesting to see whether some special assumption about the na- 
ture of P molecules would lead to such a probability of inequality (1), 
which coincides with the actual frequency of occurrence of giant chromo- 
somes. The problem undoubtedly will be complicated by considerations 
of the mechanism of multiplication of the chromosome chains themselves. 

As stated in the title, all the above is merely a collection of some sugges- 
tions for a theory. Only further research will show whether these sugges- 
tions are useful. They definitely outline a number of specific physico- 
mathematical problems which are relevant to cell biophysics. Such phe- 
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nomena as elongation of freely dividing cells, which are so well described 
quantitatively by the diffusion drag force theory, will have to be reinter- 
preted. Some conclusions of the old theory may be carried over directly. 
Equation (5) thus gives us essentially the size of a dividing nucleus, since 
complete orientation is a necessary condition for division. Though differ- 
ent from the corresponding expression obtained in the diffusion drag 
force theory, the similarity is obvious. 

If lactic acid is one of the principal sources of hydrogen ions in the cell, 
and if the latter is the cause of orientation of the P molecules, then the 
conclusions of the diffusion drag force theory about the role of glycolysis 
in cell division remain fundamentally unchanged. 

In developing a theory along the above outline it will be necessary to 
take into account the recent cytochemical studies (Caspersson, 1950; 
Leuchtenberger, 1949; Leuchtenberger and Schrader, 1951; Schrader and 
Leuchtenberger, 1951) which have important implications in regard to 
possible mechanisms of cell division. 

This work is part of the research program carried out under United 
States Public Health Service grant G-3312 to the University of Chicago. 
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The problem of social learning has been studied by psychologists and sociologists in 
the laboratory and the real world. We first consider the special, somewhat timely, 
problem of teaching a group of men to be good paratroopers. Some mathematical 
models are developed for a simple social Jearning situation, and necessary and suf- 
ficient conditions are then given for obtaining the “best” learning results or ‘‘most”’ 
learning for these models. The results are then generalized to cover other social 
learning situations; for example, classroom teaching, possibly. 


Let us first consider the group learning situation. If one learning situa- 
tion A, is “better” than the next A», we would expect the chance that all 
the persons will learn to be no less in situation A; than the chance that all 
the persons will learn in situation A;. We might also expect the chance 
that no more than one (0 or 1) person will fail to learn will be no less in 
situation A; than in A», and the chance that no more than two (0, 1, or 2) 
persons will fail to learn will be no less in situation A; than in Ap, etc. 

Now let us consider a group of men whom we want to transform into 
jumpers. By psychological and psychiatric diagnoses we can determine 
whether or not person W will be more likely to jump than person V (that 
is, in the paratrooper terminology, which is less of a quitter). We shall 

-assume that there is a certain initial chance p, (proneness to jumping) 
that person W will jump. We shall see later that it will not be required 
to estimate this probability in practice but only to say how many persons 
in the group have better chances of jumping than W. We shall say person 
W is better with respect to action J than person Z if, and only if, pu > pz. 

We expect that the event of seeing others jump may have some effect 
on a person’s proneness to jump. Let s;(p) denote the chance that a person, 
whose initial proneness was /, will jump given that he has seen 7 people 

* This paper was prepared in connection with research sponsored by the Office of Naval 
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jump. It seems reasonable to postulate that s;(p) is a nondecreasing func- 
tion of p (increasing for some values of p). The quantity 2,(p) = 5;(p)/ 
s;1(p) is a measure of the relative effect on a person, whose initial prone- 
ness was /, of seeing the jth person jump. 

Let us consider the hook-up system and corresponding model where, 
when a man who is hooked up in the airplane refuses to jump, all those 
men hooked up in line after him are unable to leave the plane. One learn- 
ing situation A, will be defined as “‘better”’ than another A, if the prob- 
ability that no more than c persons will fail to jump in situation A, is no 
less than in situation Ae, for c = 0, 1, 2,.... We might raise the ques- 
tion as to which is the best (optimal) way of lining up the men in the 
plane. What is the optimal arrangement of the men in line? 

Theorem 1: A. If there exists an optimal arrangement for the model 
under consideration it can be no better than that obtained by placing the 
best man in the first position, the next best man second, the third best 
man third, etc. 

B. The arrangement of men from best to worst is the optimal learning 
situation if and only if 7,(p) is a nonincreasing function of p. 

Proof: The chance that at least ¢ will jump in our model equals the 
chance that the first c will jump. This is, in fact, equal to 


Pr{ 72 6} = So (po) 51 (p1) S2 (pa)... Sena (Pea) ; 
where ; is the initial proneness of the person who is hooked up in the 
1+ 1 place in line. If p; > p; for some 1 > c — 1 andj <c —1, then, 
by the monotonic nature of s,(p), the total probability Pr{7 > c} will not 
be decreased if the person with the initial proneness p; changes position 
with the person with proneness p;. Hence, we have proved part A of the 
theorem. 

Now for part B. Let us consider the case A, where pth the 
--. 2 pin 2 pi 2.... It will suffice to show that this particular ar- 
rangement is better than the arrangement A» where persons i — 1 and i 
have changed positions (all others remaining in the same position) if, and 
only if, 2,(p) is a nonincreasing function of p. We have, when c > nf 


D=Pr{ 4S ¢\ Ai} —Prijscl Ay} 
= K[ss-1 (psa) $s (ps) — Sea (3) 55 (Os-3) ), 
where K is positive. Hence, the difference D > 0 if, and only if, 
[San (ps1) 53 (0s) — 53-1 (D,) 55 Oa) 1 S 0 


Si(pi) Si (Pi-1) 
Sea (Oy) See 4) 
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Since p; < pi-1, the last condition is equivalent to the condition that 
n;(p) is a nonincreasing function of p. Q.E.D. 

It should be pointed out that in the preceding theorem we had not 
dealt with the case where s,;1(p) = 0 for some values of p. This case 
causes no real difficulty and requires only slight changes in the wording of 
Theorem 1 and the results which will follow. 

Theorem 1 is sufficiently general to deal with the paratrooper hook-up 
system; i.e., where the chance of jumping is zero after there has been a 
failure in the plane. However, in general social learning situations it is 
clear that the event of seeing others fail may have some effect on a per- 
son’s initial proneness to perform an event, but clearly the effect will not 
be to reduce his chance of performance to zero. Now let s;,(p) denote the 
chance that person, whose initial proneness was #, will jump given that 
he has seen 7 people jump and & people fail to jump. It seems reasonable to 
postulate that sj.(p) > sj, n(p). 

We now have the following intuitively clear 

Lemma 1: The probability of obtaining at least d further failures given 
that there have been jumps and & failures is no greater than the probabil- 
ity of obtaining at least d further failures given that there have been 
j — 1 jumps and & + 1 failures. ' 

Proof: For d = 0, we have 

Pr{ fsd=0| 39, BY = 55, (Psy) Spar, y Pyrnts) Sie, w Pyrete) so - > 
hence 

Pr{ f<d=0|7, R}>Pr{f<d=0lj—-1, k+1}. 

Now let us consider d = 1. The lemma is true when Pr{f < d|j,k} =P 
is a nondecreasing function of Smn(p) = s. We therefore wish to show that 
5P/5s:> 0 (we have assumed differentiability of s to simplify the presenta- 
tion, but this is not a necessary assumption). For d = 1, we have 


BS 2p 944, bO742, Be AOE ko 74h, k° LRG ee os OG is 


and 
>a Snipe, ko?" CES) Sg ur mene re Sa eee 
when n=k+1 
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Hence 
EOP >0 
6S. k+1 
and 
6P 
‘sccm eal Gog EVE, pene Sim—1, eS mt, w>mt2, ko °° 


5 bit na lke ee 


where K is non-negative. We also have that 


bP _ > 


oS) 7 


BUNCE: (Spt, Esmee? isco ce Sate oe eee ee the preceding result 
for d = 0. We may now proceed to prove the lemma in this way by 
mathematical induction on d. Q.EL.D. 

In analogy to the hook-up model, we define 


5j,%(P) 
S72 EP) 


where f; , = 1 — s; x. We have 

Theorem 2: The arrangement of men from best to worst is the optimal 
arrangement if 7j.(p) and q;.(p) are nonincreasing functions of p. 

Proof: Let us consider the case A; where pp>fithae...2 
pia > pi >....It will suffice to show that this particular arrangement 
is better than the arrangement A» where personsz — 1 andz have changed 
positions (all others remaining in the same position). If 7 jumps and i — 
j — 1 failures have been observed by the person in the 7th position, the 
possible results for the ith and i + 1st person have the following prob- 
abilities for case Ai: 


fin (pb) 


ane 9: (B) = fae 


N jk (p) to 


(Opt Liege (P,_1) Saha ee seal (P,) 
yi Sco OS) Spee PD een Pe eee 
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If we now consider case Az, where the p; and p;_; are interchanged, and 
then subtract the probabilities obtained, we have 


{0},— {0}, = Snipa es) Sha sed ee 
ny ee (2;) a ee (P,_,) 


= 1M, i (P;) a oa, i= 2s) j K20 


SOCIAL LEARNING SITUATIONS Shall 


and 
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We write {0}1 — {O}2 = do(7, i -—7 — 1) and {2}. — {2}. = (7,7 — 
j — 1). We have then {1}; — {1}. = &(j,i —7 — 1) — do j,i —7 — 1). 

Now let us consider the probability Po that all will jump. There is no 
difference for this case between the hook-up model of Theorem 1 or the 
present model. Hence, it is clear that the probability that all will jump 
is no less for A; than Ay. We have, in fact, that the difference between Po 
for arrangement A; and P, for arrangement Az is APy = (0) + doi — 1, 
0) - [0, 0], where (2) is the probability that there will be & failures among 
the persons who perform before person z — 1, and [k, ¢] is the probability 
that there will be ¢ failures among the persons who are to perform after 
person 2, given that there have been a total of & failures among the per- 
sone Oole2, . tat sh, 1. 

We are now ready to consider the probability P; that one will fail to 
jump. We may differentiate three cases as follows: the failure occurring 
before person 7 — 1, after person 7, or by either person z — 1 or z. Hence 
we have 


AP, = (0) -do(i—1, 0) + [0, 1] + (0) - [d, @—1, 0) 
—d)(t—1, 0)]-[1, 0] + (1) +d) @—2,1)-[1, 0], 
and 
AP) + AP; = (0) +d) (i—1, 0) -{ [0, 0] + [0, 1] — [1,0] } 
+ (0) -d2(i—1, 0) + [1, 0] + (1) +d) @—2, 1)-[1, 0]. 


From Lemma 1 we have that [0, 0] > [1, 0], and, therefore, AP) + AP; > 
0. 
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Now let us consider the probability P. that two will fail to jump. We 
have 


AP: = (0) +d) (2—1, 0) - [0, 2] + (0) + [d2 @— 1, 0) 
—d)@—1, 0))<11, 1)-+ (1) do e251) ae 
— (0) +d, @= 1,0) 12) 0) CL) hs 2 


—= "dg (tae 1)]-[2, 0] + (2) od i —3, 2) 1 (2; 0] 
and 


AP) +APi+AP2 = (0) +d) (i—1, 0) + { [0, 0] + [0, 1] + [0, 2] 
~ [1, 0] S-[1,. 1] } + (0) 2d, G@— TO it Ole | et eal 
+ (1) +d) @—2, 1)-{ (1, 0]-+ ia a2 Oe 
+ (1) +d2(i— 2, 1) + [2, 0] + (2) «dy (G—3, 2) + [2, 0]. 


We may now proceed to prove the theorem in this way by mathematical 
induction on the number of failures. Q.E.D. 
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In this note the explicit solution is given to an equation, suggested by C. N. Hin- 
shelwood (1946), describing the growth of a bacterial population under the assump- 
tion that toxic products are a limiting factor. The behavior of the culture as a function 
of time and the parameters (initial number, rate of growth, and rate of production of 
toxic substance) is discussed. 


C. N. Hinshelwood (1946) in his studies on bacterial growth processes 
has considered the relation between the number of organisms in the cul- 
ture and the concentration of some toxic metabolic product. Hinshel- 
wood’s equation, which states that the number of organisms will increase 
at a rate proportional to the number at any time and decrease at a rate 
(per organism) proportional to the concentration of the toxic substance, is 
as follows: 

Ot = bn (1— ac); (1) 
where is the number of organisms in the culture, c is the concentration 
of the toxic metabolic products, and & and @ are positive constants. If the 


toxic substance is formed at a rate 7 per organism, we have 


dc _ oe 
aha 
or } 
= nat. (2) 
(a rf 


We can now write (1) as 


dn : 
at = kn (1-ar f nat). (3) 


The purpose of this note is to obtain an explicit solution for m(¢) and dis- 
cuss its behavior as a function of the parameters. 


* Public Health Service Research Fellow of the National Cancer Institute, National Insti- 
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To solve the integrodifferential equation (3) we divide both sides by ” 
and make the substitution y = In 2 or = e”. We have now for (3) 


biid = —ake. (4) 


We have put a = ar in (4). Multiplying (4) by 2dy/dt and integrating, 
we have 
CLE 
(3 a eae 


where C = k? + 2akmo, since n = mp and dy/dt = k when ¢ = 0. Now 


D = + VET Dak (iy &)- (5) 


Expressing (5) in terms of 7 we have 


ed wi ai pe 
— a VE + 2ak (iy — 2). (6) 


Thus for the range of m where dn/dt is positive 
2 dn 
lal Reema ; Mg 


Performing the integration we have 


_1, [(VAt=2akn— A) (k+A) 
mre versrsrar ae eee | (8) 
or, solving for n(t), 
G At 


In the above we have put, 


A= VRP+2aknm, 
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and 


We begin the study of the behavior of n(#) by differentiating (9) to 


obtain 


on AC e4* (1 —Be4t)-24 2 ABC e24* (1 —Bet)-3, 


Setting the derivative equal to zero and solving for ¢, we obtain 


1 1 1 k+A 
a —_ — = —_— 
t zn ( a) =n ( =) (10) 


Differentiating again and introducing ¢* for t, we find that the second 
derivative is less than zero. Hence n/(t) has its maximum value when f 


n(t) 


Nm - 


FicureE 1 


assumes the value given by (10). By introducing ¢* in (9), we find the 
maximum number of organisms to be 


Ny = Mot. (11) 


In (10) we have the time required for the population to reach its maximum 
value as a function of the initial number of organisms, and in (11) we give 
the maximum number of organisms as a function of the initial number. 

Beyond the maximum point d/dt is negative, and a completely analo- 
gous calculation shows the growth curve to be symmetrical with respect 
to the line ¢ = i*. The complete growth curve is shown in Figure 1. 
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It is of interest to point out that (9) can be considered as a modification 
of the well-known logistic equation which occurs in many population 
studies. The distinguishing feature of our equation is that the denominator 
is squared, while in the logistic equation it is not. Moreover, in the logistic 
equation B = —C. 

The author wishes to thank Dr. Anatol Rapoport for suggesting the 
problem and for helpful discussions. 
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A neural mechanism is studied in which three primary receptors interact in such a 
manner that only a small region is in activity for a given colored stimulus. The 
similarity to the color triangle is noted. The introduction of an additional primary 
receptor may or may not change the dimensionality of sensation space in this 
mechanism. Problems of discrimination of color, absolute judgment of saturation, 
and hue are discussed. 


It is the purpose of the present paper to develop some aspects of a 
neural mechanism which were introduced to represent some qualitative 
features of color vision (Householder and Landahl, 1944). It was assumed 
that there are three types of primary receptors, responding differently to 
various wave lengths A of light, in such a manner as to cover the visible 
range. Let p;(A) be a measure of the response of the ith receptor to a unit 
intensity of wave length d. For a given light stimulus with an energy 
spectrum S(A) let there result a stimulus of intensity 


Si = JS (d) p1 (A) dd 


in the first type of receptors. The stimulus intensities S. and S3 of the 
other type receptors are similarly defined. To a given S(A) there is there- 
fore a corresponding excitation pattern which we may denote by (Si, So, 
S3), or, more simply, by (S;). 

Let each receptor of type i be connected to neural elements lying in a 
plane C defined by three points C:, C2, C3. We consider the region C to be 
a plane for convenience only. Let the number of neural elements in C be 
large enough so that we may speak of their density per unit area, and, 
similarly, let the number of connections from the type 7 receptors be large. 
We now define p;(x, ) as the density of connections from type 7 receptors 
to some point (x, y) in C. For convenience let Ci, C2, and Cs be equi- 
distant from each other, the center of the triangle C,C2C being the origin 
(0, 0). Let the x-axis be taken in the direction of Ci, the units being chosen 
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such that C, is at (1, 0) and thus C2 is at (—3, cos 30°) and C3 at (—3, 
—cos 30°). 

We shall now suppose that each of the density functions pix, ¥) de- 
pends only upon the distance d; from the point (x, y) to the corresponding 
Cy, i.e., pi(x, v) = o:(d,), and that o(d) is a monotonic decreasing function. 
It might be supposed that o(d) is a normal error distribution due to some 
random process in development and that since we shall only be interested 
in values of d which are small, so that if oo? and y; are constants, o;(d;) is 
given by 

o: (di) = (1 — vedi) 00: « (1) 


It can be shown that the replacement of d? by |d;| leads to similar results. 

The intensity of excitation EZ; (Rashevsky, 1948) is a function of the 
intensity S;. We shall consider two cases, one in which the relation is linear 
and the other in which it is logarithmic in the external intensity S;. Thus 
in the first case we have, if /; is the threshold, 


E; = a; (Si— hi). (2) 


If we consider only the steady state, then the density of the excitation 
factor e; at any point (a, y) due to EZ; is given by E;o; so that we have 


e;, =Ljo; = dog; (S;— hi) (1 — yd) (3) 
and, hence, if we define 
Si = ajo: (Si— hi) =Ejos , (4) 


the total excitation density (x, y) is given by 


e(4, ¥) =atete =S1+534+S3— ySidi — yoS2d3 — y2S3d3. (5) 


If the ’s are about equal and not too large, e(x, y) is a function having a 
maximum value within the triangle C:C,C;. Suppose now that the region 
C is connected with another region C’ by excitatory pathways, in such a 
way that there is a point to point correspondence between the two regions, 
preserving the geometry. Let there also be numerous cross-inhibitory 
paths from each connection in C to every pathway originating in C’, 
which thus tend to restrict the excitation in the region C’ to those ele- 
ments in the neighborhood of the point P’ of C’ corresponding to the P 
of C for which e(«, y) is a maximum (Rashevsky, 1948, p. 377). Thus for a 
particular S(\) only a certain region in C’ will be excited and the position 
of this region will depend upon S(a), i.e., the relative amounts of energies 


COLOR VISION 319 


of the different wave lengths of the incoming light. Any S (A) giving the 

same values for S;, S:, and .S3 results in the same region being excited. 
We may solve for the coordinates of the center of the region as a func- 

tion of S;, S, and S3 by setting d¢/dx = d/dy = 0, and thus obtain 


Se 2S me v2.59 ; 353 (6) 
2 (7,5! + 7,50 1,5) ° 


/3 (-y25'3 = 353) (7) 
2AVS7 ey Sock VIO) 


es 


If we set X* = x*, V* = y*V/} — x*/3, Z* = —x*/2 — V/3y*/2, X, 
Y, and Z become the coordinates of a point in the coordinate system OC;, 
OC,, OC3. Note that X + Y + Z = Oand Y and Z can be obtained from 
equation (6) by cyclic permutation of the subscripts. The point C; can 
be represented by (1, —3, —}), C2 by-(—4, 1, —4), and C3 by (—3, 
—%, 1). If now 1S; > y2S; and 353, then X = 1, Y = Z = —3, 50 that 
only those neurons in the region about C{, corresponding to Cj, are active. 
From symmetry it can be seen that if y253>> Si, 7353, only neurons 
about C3 are active. On the other hand, if y,S{ ~ 253 >> 353, the ac- 
tive region in C’ lies near the line between C; and Cj, while if 5; = 
2S = 7353, the active region of C’ is at the center of the triangle formed 
by Ci, C2, and C3. Note that for any value 8, if all 6S; and S; are far above 
threshold, the stimulus with an excitation pattern (051, 0S), 653) results 
in an excitation with a maximum in almost exactly the same position in 
C’ as does (Sj, So, S3). But if any S; or 0S; is near /;, a general change in 
intensity by a factor 6 changes the position in C’ by an amount which 
depends upon the relative values of the /;. Furthermore, if S$; >> h; for 
eachz, and if (S;) = (Si, 52,53), the center of the region which corresponds 
to the stimulus (S; + 0/aio0,;7;) is always nearer to the origin than the 
point which corresponds to (S;), the point corresponding to (S;+0/aicorv:) 
being at the origin if 6 is large enough. 

If Sj = Sz; = 53, and S; > hz, the intensities 4, @, and €; at C1, C2, and 
C3 respectively are equal. Also if yi = Y2 = 7s, these physoloelealy 
equal values for the e; result in an excitation about the origin of CH the 
distributions o; are the result of the same random process we might expect 
this would be so and y; = Y;; also oo: = 0;. In most of the following we 
shall make this an assumption and thus drop the subscript ; of oo; and 7. 

From the above discussion it can be seen that we have mapped the 
familiar color triangle (cf. Wright and Martin, 1946) on to the neural 
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region C’. Energies impinging upon three distinct receptors result in a 
response somewhere in a plane region, the position depending upon the 
relative amounts of energy in the three components. The primaries them- 
selves have no special properties, and a mixture of any two does not pos- 
sess any special properties of either. 

A change in S(A), which causes the corresponding point P’ (in C’) to 
move radially, corresponds to no change in hue. This type of constancy 
can be taken into account by the following addition to the mechanism. 
Let there be neuronal elements lying along some closed curve Cy con- 
taining the triangle C,C,C3. Let each neuron of Cy have dendritic connec- 
tions with the axon endings of those neurons of C’ in such a way that con- 
nections are made only radially toward the origin O’ of C’ and in such a 
manner that the probability of a connection decreases to zero as the 
origin is approached. Then the neurons in the same neighborhood of Cy 
are excited when there is no change in hue, the excitation being greater, 
however, for a more saturated color. The assumption that the connections 
of the neurons in a small neighborhood of some point of Cy be made al- 
most exactly radially is not at all necessary. There can be a very con- 
siderable lateral spread and even a spread past the origin if Cy transmits 
through to another Cz only in the region of maximum excitation. Such a 
radial structure as that of the center Cy requires only a simple tropism for 
its development. 

In a similar manner, a center Cg, in which the neurons in a linear array 
each receive impulses from different annular regions, would result in a 
response which would depend almost entirely upon the degree of satura- 
tion. Such a structure could arise as a result of the various elements of 
Cs being attracted, during development, to various absolute concentra- 
tions of some substance produced at the origin of C’ which then diffuses 
outward. Thus C's is a much less probable structure than Cy. It may be 
noted that absolute judgments of saturation are relatively poor compared 
with judgments of hue. 

We shall next determine the size and shape of a region in C’ which be- 
comes active for a given S(d). Two colored stimuli producing excitation in 
regions which just fail to touch should correspond to two stimuli which 
are just clearly distinguished. Thus the radius R of the active region in 
C’ would be a suitable measure of the psychological distance between two 
colored stimuli. The radius R is small where the discrimination is fine. 
Suppose then that transmission from C to C’ occurs for values of ¢ within 
an interval émax > € > €max — Ae and the range Ae is nearly independent 
of the total intensity of illumination. Then, since the radius of curvature 
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of e(x, y) can be shown to be constant, we find the radius R of the active 
region of C’ to be given by 


a (8) 


1 (SESS? 


so that accuracy of discrimination increases slowly with illumination. For 
not too large distances, points on the triangle Cj, C2, C3 which are equally 
distant are distinguished with equal difficulty. 

If, however, one uses a secondary set of “primaries,” the circles implied 
by (8) will become ellipses. Furthermore, if one plots a circle in the C’ 
plane, the result is an ellipse in the plane f;f; where f; is defined by 


we Qiyi (S; a h;) O04 
“" Dasv3 (S3— hy) 095° 


If S;>> h; for each 7, the quantity f; becomes the fraction of energy going 
to 7 in arbitrary units. In this case a pair of f’s may be chosen as coordi- 
nates, the results of any plot depending upon which pair is chosen. A plot 
of this kind is conventionally used, and it is on such a plot that the data 
yield ellipses (cf. Brown, 1951). If one considers the more likely possi- 
bility that the points Ci, C2, and C3 are not exactly equidistant, the result 
is further complicated. 

In the above case we have assumed F; to be linear in S;. This is gen- 
erally not the case for visual phenomena. It is more often assumed, as a 
result of various observations, that the relation is logarithmic. We shall 
now consider the effect of this complication. Instead of equation (2), we 
will now have, if 6; is a constant, 

Si 

E; = 8B; log 7 (9) 
If no further changes are made, we find that in this case the active region © 
of C’ changes its position slowly with change of intensity, the proportions 
of the S; being unchanged. Generally the change is such that at higher 
intensities the position of the active region is nearer to the origin. The 
change appears to be too rapid to correspond to the actual situation (cf.. 
Brown, 1951). We shall therefore consider another mechanism. 

Let each primary receptor produce an amount j of inhibitory factor 
(instead of excitatory) at a connection A; along the pathway to C. Let 
there also be excitatory elements which react to the total physiological 
intensity (1/1 + S2/d + S3/hs) and which act at each of the connec- 
tions A; along the pathways between the receptors and C. Then the 
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amount of excitation e; at A; will be the quantity (’ log (Si/d + So/In + 
S3/hs)/h’ — B log Si/i. Then e will be given by the following expression, 
instead of by (3), 


Si, Se, Sa 
Ta a hs Si 
(x, 9) - SS ou nnn pin 3 *— Blog 5*|, (10) 


the d,’s being the same distances as in equation (5) and thus being func- 
tions of x and y. If we simplify matters by assuming h’ = 1, B; = 8’, 
ViOoi = Yj00;, and introducing the fractions fi = Si(Si/In + S2/I2 + 
S3/h3)h:, so that fi + f2 + fs = 1, then we find the position of the maxi- 
mum for any Sj, S2, Sz to be given by 


Joe 22/3 
og nel 
ae (11) 


2 log 


with similar expressions for Y and Z obtained by permuting the subscripts. 
In this case, the position does not depend upon the level of the intensity 
but depends only upon the relative intensities of the components. If 
fix 1 or S;> So, Ss, then X = —§ and Y = Z = — 3; while if fi = 
O(S1 >> Se, 53), X ~ 1 and Y = Z = —}3. Thus the representation is in- 
verted with respect to the linear case. The radius of curvature in this case, 
instead of being given by (8), is a constant, independent of «, y, and S;. 
To study the differences between the two cases corresponding to equa- 
tions (6) and (11), we may (since in each case the radius of curvature is 
independent of x, y) compare the effects of changing the relative value of 
one of the primaries, e.g., fi. Consider a change such that f. = fs; then 
form the expression Af;/Ax™. This expression shows how great a change in 
fi is required to produce a given change in x subject to the condition that 
fx = fs. For the first case (linear, A, Fig. 1) this ratio is found to be equal 
o 3 if f; is defined by 
a; (S;— h;) 
Za; (Sj— =H) 


fi= 


In the second case, corresponding to equation (11), we find (B, Fig. 1) 


Afi _ 1 fr —f) [log fi+2 log (1 — fi) — log 4]? 
Ax* 3 (1— fi) log 2— fi log fr— (1 — fi) log (1 — fy) * 


(12) 
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In a similar way we can obtain the change Ay* in y* when fy and fs = fo 
are changed to f, + Afe and f; — Afy. We then find Af,/Ay* = 1/V3 in 
the linear case (C, Fig. 1) while in the second case (D, Fig. 1) we obtain 


Af, _ (=f) 
Apt = TAF log fa 2 log (1 fa) —2 log 2]. 13) 


The signs are reversed in the curves of Figure 1 for the case correspond- 
ing to equation (11), since, as was pointed out above, the representation 


2 4 6 8 1.0 
ene 


Ficure 1 


in this case is inverted with respect to the first case. In the first case the 
change in # and y directions are nearly alike and are constant. In the sec- 
ond case a primary and its contrasting color are dissimilar. A small addi- 
tion of a contrasting color A(1 — fi) toa large amount of a primary (f=1) 
is relatively ineffective in producing a noticeable change in response (cf. 
solid curve for f ~ 1 in Fig. 1). On the other hand, only a small amount 
of a primary Af; added to a contrasting color (f2 = f3 > fi) is necessary to 
produce a change in the response (solid curve for fi ~ 0). 

We shall now study some effects which might arise if more than three 
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primaries are involved. We shall only study the simpler linear case as this 
illustrates the effect of the introduction of additional receptors. 

Let there be four kinds of primary receptors which send out pathways 
into the region C in a manner similar to that considered in the first case 
above. First suppose these pathways spread out from four points Ci, . . . 
C, which are the corners of the square CiC2C3C4. Let the distribution 
functions of the densities of connections o; be given as above by equation 
(1). We can then find the point of maximum activity as above, the result 

oO 
being ere 


ree meE REC TE, aay 


with expressions for Y, Z, and W obtained by permuting the subscripts. 
However, since the points lie in a plane, mixtures of only three stimuli can 
be made equivalent to any compound stimulus which lies within the tri- 
angle formed by the three stimuli. Thus from the data on mixture of 
colors alone, the representation by three ‘“‘primaries’”’ would be adequate. 

If this case were ‘‘normal,’’ then an abnormal case arising as a result of 
the absence of one of the primaries (or its failure to act at C or its acting 
uniformly over C) would appear to be almost normal so far as color mix- 
ture is concerned. However, in the case of three primaries, a complete loss 
of one would radically change the response to color mixtures. In this case 
all colors would be arranged along a line. Note that with a partial mal- 
function of one of the primaries it is possible for the representation to be 
well spread out over an area or to be essentially a linear arrangement or 
any intermediate situation. 

Now if there are four primaries, but the structure is such that the points 
Ci, ... C4 are not coplanar but, for example, form a regular tetrahedron, 
then we find 
= (351s ee) 


A 35 SESS: 


with similar expressions for Y, Z, and W. In this case it is evident that a 
color match would require the mixture of four colors in general. 

If the functions ,(\) are such that, for a moderately high intensity, 
Si > hi, t = 1, 2, 3 regardless of S(d), then there is a band of value near 
the border of the region C’ which would not normally become excited. 
The points in this border region would correspond to highly saturated 
colors. For example, if a particular S(A), with given total intensity, results 
in S;/S; having a maximum for S; = S3 > = h3, then x < (S; — 
S2)/(S1 + 252) = x* < 1 and no color can produce excitation at a point 
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«x > «*. However, if the receptors of type 2 and 3 are partially or com- 
pletely inhibited in any manner, then the result is the activity in the 
region about « > «*. This inhibition might arise as the result of accom- 
modation to a color S$; ~ S3>> S; or due to inhibition from an adjacent 
region in which S; = S;>> S,, provided that the cross-inhibitory effect 
results in a uniform reduction in all S;’s or S; in region J inhibits only S; 
in region JJ, and conversely. In any such situation one might expect to 
obtain colors appearing to be more saturated than would otherwise be 
obtained. 

It may be noted that in the case of three primary receptors, the inten- 
sity can be included so that the result is a color pyramid (Householder and 
Landahl, 1945). On the other hand, the intensity may be given independ- 
ently. The results of experiments on color discrimination will be different 
in these two cases, the representation of intensity having similar geometric 
properties to color in the first case, while having different properties in the 
latter. A careful comparison with experimental data may show one ar- 
rangement to be a better model than the other. Available data (Brown, 
1951) do not appear to be adequate to shed light on this point. 

At high enough intensities, the functions E; must physiologically be 
_ bounded so in either case, when the approximations (2) or (9) break 
down, the active region in C’ will, regardless of the stimulus, move in to- 
ward some point located not far from the center. 

In most of the above discussion we have considered three groups of 
receptors, those in each group being considered to be identical. It can be 
‘shown that to account for most of the observed data on reaction to color 
it is only necessary that there be three groups of receptors with the proper 
geometrical connections, but each having average functions (A) suf- 
ficiently distinct. 

The author is indebted to Dr. George Karreman and Mr. C.S. Patlak 
for reading and discussing the manuscript. 
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If, along the length of a blood vessel, there is a discontinuous change of diameter, 
the pressure waves will be reflected at the discontinuity. The transmission and re- 
flection coefficients of the waves are calculated in terms of the ratio of the diameters 
and the ratio of the propagation velocities of the waves in the parts of the blood vessel. 
We obtain also the reflection coefficient in terms of the ratios of the diameters, of the 
elasticity moduli, and of the wall thicknesses. 

The possible clinical importance of the results is pointed out for cases of coarcta- 
tion, arteriosclerosis, and other pathological conditions which might effect local 
changes in the structure of the blood vessel wall. 


Whereas the propagation of pressure waves in a fluid within an elastic 
tube has already been studied (see Boulanger, 1913; Witzig, 1914), their 
reflection has not been studied so extensively. W. F. Hamilton and P. Dow 
(1938) claim to have obtained experimental evidence for the existence of 
standing waves of the blood pressure in the aorta. In another work Hamil- 
ton (1944) has shown that general vasodilatation due to the injection of 
acetylcholine caused the reflected waves to become weaker. With larger 
doses of acetylcholine the standing waves were eliminated from all ar- 
teries. That the drug had not affected a change in the arterial wall which 
prevented the transmission of the reflected waves was shown by the fact 
that these waves reappeared when the leg was grasped tightly at the knee 
producing a vasoconstriction. This, together with the observation that 
they also could be brought back by intra-arterial injections of epinephrine 
in small doses producing vasoconstriction of the dilated arterioles, strongly 
suggests that the change in radius of the arterioles plays an important role 
in the production of the reflected waves. This is shown in the purest form 
by I. G. Porjé (1946), who concludes that in cases of patients with a 
coarctation (isthmus stenosis) some of the waves are reflected more power- 
fully than in normal individuals. Porjé also emphasizes the importance ofa 
theoretical study of reflections of pressure waves in the arterial system. 
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The lack of a theory which includes the study of the reflection relations of 
the wave has been, in his opinion, the great obstacle in the study of the 
pulse wave. A similar opinion was expressed some years ago by A. Apéria 
(1940), whom Porjé quotes. The presence or absence of the reflection of 
pressure waves gives some indication of the cause of a low blood pressure 
and the degree of vasodilatation in that condition (Hamilton, 1944). If the 
blood pressure falls due to a hemorrhage, the arterioles constrict and there 
will be strongly reflected waves. However, if the blood pressure is low be- 
cause of dilatation of the peripheral arterioles, reflected waves are absent. 
The evidence given above suggests that a change in the diameter of the 


A B 


Ficure 1 


blood vessels probably plays an important role as far as the reflection of 
pressure waves in that vessel is concerned. Therefore, we will investigate 
theoretically the simplest possible model—directly suggested by the coarc- 
tation studies mentioned above—that of a discontinuous change in the 
diameter of the vessel (Fig. 1). 

As at any discontinuity a pressure wave traveling from A to B is par- 
tially reflected at the sudden change BC and partially transmitted. To de- 
rive expressions for the reflection and transmission coefficients we have to 
take into account the two conditions (Rayleigh, 1896, p. 63): 


Pa=Pc, (1) 
the continuity of the pressure, and 
tRvp=TPrv¢, (2) 


the continuity of the mass flow, v representing the velocity. 

Actually, instead of (1), the Bernouilli equation should be used, but, by 
taking into consideration that the velocity heads are negligible compared 
with the pressure heads, (1) results. 
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Any pressure wave traveling without deformation with velocity ¢ from 


A to B can be described by 
x 
p=F(t—2), (3) 


where ¢ represents the time, x the abscissa measured from an arbitrary 
origin, as, e.g., A, and F denotes a symbol for a function of the argument 
t —x/c. It can be easily seen that this is the case. Consider a pressure 
disturbance, or, what amounts to the same thing in elastic tubes, a wall 
disturbance traveling with velocity c to the right. Suppose at ¢ = 0 it is 
at x = 0, so that 


i— 


<= 0. 
Cc 


Corresponding to (3) the pressure at x = 0, when ¢ = 0, is then p = F(0). 
At the time # the disturbance traveling with velocity c is then at x = ct. 
Therefore, that same pressure F(0) is at ¢ at x = ct, and then, because 


we have, according to (3), = F(0) which proves our statement above. 
A wave traveling in the direction of the negative x-axis, that is, from 
B to A, can be correspondingly represented either by 


x 
eee): 
because the velocity is then —c instead of c, or by 
if, 
= = 4 
p=s(t+4). (4) 


Let us assume now that a pressure wave represented by (3) is partially 
reflected at the discontinuity BC producing a reflected wave represented 
by (4) and a transmitted wave which, similar to (3), may be represented 


e p=s('-4), (5) 


the velocity c’ in the tube CD being usually different from the velocity C 
in AB. The total pressure in AB is then the sum of the pressures given in 


(3) and (4): per (eS (6) 


whereas the pressure in CD is represented by (5). 
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For convenience we will take as our origin in the remaining part of this 
paper the discontinuity BC. Then we have 


pe=Paz,,~=FuM+f®, (7) 


whereas 
Poe Pcie, ee (8) 


Introducing (7) and (8) into (1) leads to our first relation between the 
incident, reflected, and transmitted waves at x = 0: 


Pit) Oa ree (9) 


We will now consider the velocity which is connected with the pressure 
wave (3). This relation could be obtained directly by integration of the 
Euler equation, using (3), but we will do it ina more elementary way. To a 
small element of fluid in a tube with radius 7 and length (Ax) we will apply 
Newton’s law 
K=ma, (10) 

in which K represents a force which gives to a mass m an acceleration a. 
The volume of that element is rr?Ax. The mass m is therefore prr?Ax, 
p being the density. 

If the velocity is represented by v(¢), the acceleration a is, on the aver- 
age, 
t+ At) — 2) 

At i 

The force K on the element is the difference between the force due to the 
pressure on the left-hand side, 


p (x) ‘ila ) 
and that on the right-hand side, 


p(x +Ax) rr?. 
Therefore 


K={p(«%) —p(a+Ax) }rr?. 


Introducing these expressions for the mass, acceleration, and force into 
(10), we obtain 

{p (x) =p (x ++ Ax) lar? = prr’Ax eae) 3 
or, simplifying, 


p(x) — p(x+Ax) 3 v (t+At) — 9 (t) 


which is an approximate form of the Euler equation mentioned above. 
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If p is given by (3) we have, introducing (3) into (11) and rearranging, 


a ee 
Ub At) 10.) = z 
Taso i) - Ag . (12) 
Alt 
The denominator Ax/At is the velocity c, 
Ax | 
=o: (13) 
from which we find 
Ax 
— = At. (14) 


Introducing (13) and (14) into (12) we obtain after some new arrange- 


ments: 
v(t-+At) — v(Z) == \F(t-*)-F(s-ai—*)h, 


or 


v (t+ At) — 0(2) -—F(} (t— At) ~=t+at| 


-[-afe—a0 2H] 


From (15) we see, particularly if we write ¢ instead of (¢ — Af) in the 
second member, that 


(15) 


(2) -—F (1-4), (16) 


which is the relation we have been looking for. Correspondingly, the 
velocity associated with (4) is given by 


1 x 
0) = ¢ (1+) (17) 
and that associated with (5) by 
1 x 
v (2) == 8 (1-4). | (18) 


The total velocity in AB associated with the pressure given in (6) is then 


sus=aoF (1-2)——- (444), (19) 


whereas the velocity in CD is given by (18). Equations (6) and (19) are 
the basic equations of the theory of the water-hammer phenomenon 


(Allievi, 1903). 
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Then we have 


i= babe =F O -—f (t) (20) 
and 
Yo = PeDsay =F Gl g(t). (21) 
Introducing (20) and (21) into (2), we obtain: 
{Py —pbaarts eo f 
or, simplifying, 
FO -1H=55e. (22) 


The equations (9) and (22) give us two linear equations for the two un- 
known ratios f(t)/F(é) and g()/F(@). 
Solving (9) and (22) for these unknowns, we find 


2 


ieee 
POO ee 
F (i) 2 ¢ 
TRG 


or, after simplification, 
f(D SUReC-= 6 
FQ) Rete” 


(23) 


Similarly, 
gt) = 2R? c’ 
F(t) Réc’+re 


(24) 


It may look paradoxical that the left-hand sides of (23) and (24) seem to 
depend on #, whereas the right-hand sides are independent of t; the veloci- 
ties depend in general on the geometrical and elastic properties of the tube, 
as we will see, but not on ¢t. However, in the case where 


P(t-2)= 4 sinw (1-*), (25) 
s(1+=)=B sin} @ (1+=)+ot, (26) 


g(1-2)=Dsinfo(1—*) +4, (27) 


in which # and y represent the phase differences of the incident, reflected, 
and transmitted waves, we have, by introducing (25), (26), and (27) into 
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(9) and (22), after substituting « = 0 in the last two equations, and de- 
noting, for abbreviation, (r?/R?) (c/c’) =k: 

A sin wi+ B sin (wt+ ¢) = D sin (wt +) (28) 
and 

A sin wt —B sin (wt+ ¢) =kD sin (wity). (29) 


Because (1), (2), (9), and (22), and therefore (28) and (29), hold for all val- 
ues of #, the coefficients of sin wt and those of cos wt in both sides of each 
equation must be equal, which leads to: 


a --B cos ¢ = D cosy; (30) 
Bsin d= Dsiny , C3 1) 
A—Bcos¢=kDcosy, (32) 
—Bsing=kDsiny. (33) 


These four equations must be solved to find the four unknowns B/A, 
D/A, $, and w. Because is in general not equal to —1, (31) and (33) can 
only be simultaneously satisfied if 


sin @¢ = 0 (3 4) 
and 
sny=0. (3:5)) 
Equation (34) is satisfied if 
og=mr, m=0,+1,... (36) 
and (35) if 
y=nr, (ee eee bes ot ear (37) 


corresponding to a reflection and transmission in or out of phase. But for 
these values of ¢ and y the ratios 


f(t) _Bsin(wi+¢) 
Pye Asin of 


and 
g(t) _ Dsin(wt+y) 


F@ «Asin ot 


are independent of ¢ and, except for the sign, equal to B /A and D/A, the 
ratio of the amplitudes. This is the reason that we will call those ratios 
the reflection and the transmission coefficient, respectively. The reflection 
coefficient is then given by (23) in terms of the radii and the velocities, the 
transmission coefficient similarly by (24). To get those coefficients in 
terms of only the geometrical and elastic properties we have to express the 
velocities in terms of those quantities. Here we could use the results of the 
more refined treatment, given in the Appendix, which takes into account 
the viscous properties of the fluid and the elastic properties of the wall, 
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both in some detail. We will also give here an elementary treatment which, 
though crude, exhibits all the features essential for our present purposes. 
However, before doing that we would like to point out an extension of 
formulas (23) and (24) for the case in which there is a branching of a vessel 
(Fig. 2). 

The result for the reflection coefficient, in which we are here mostly 
interested, is: 


jet Bade 
EOes 2 Grehe Co (38) 
PF t) r: Cc re Cc 
3! 2 
Ve ae es 


The generalization of this formula to the case of any number of branchings 
is obvious. 


-_ > 
| 
a 


< 


| 


FIGURE 2 


To derive the relation between the velocity c of wave propagation and 
the geometrical and elastic properties of the system we proceed as follows. 
First we derive the relation between the pressure # in an elastic tube and 
the tension o in the wall (Fig. 3). The force on a piece of the tube with a 
length LZ (Fig. 3) due to the pressure on the inner wall of that piece is, 
according to a well-known theorem of physics, 


pXarea ABCD = pX2rL. (39) 
The force on that same piece due to the tension in the wall is (Fig. 3): 
o Xarea(ADIH+BCGF) =o X (6L+ 6L) =2c8L. (40) 
In equilibrium the last two forces are equal, therefore 


2rLp=26Le 
or 


¥ 
ae a (41) 
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From (41) we see that if the pressure is changed from the value p() to a 
value p(¢ + At) the tension o is changed from a(t), given by (41), if we 
take p = p(é) to a value o(t + Ad), which is correspondingly given by* 


o (t+ At) =< p (t-+ Ad). 
The change in tension is, therefore, 


o (¢+At) —o (2) =={p(t+Ad) =p) } (42) 


FIGURE 3 


Assuming Hooke’s law to hold, we have 


POAD — 7 @) =o GAD. ce @) eu (43) 


in which E£ is the elasticity modulus. 
Introducing (42) into (43), we obtain 


tba) —r() ={p +a) —p OIE. (44) 


The corresponding increase in volume of an element with length Az is 
given by, using (44), 
Zar 


Qar{r(t+At) — 7 (t) jAx= cE Ax{p(t+At) — p(t) }. (45) 


* Here we neglect changes in r and 6 as is usually done in the water-hammer theory (cf. 
Jaeger, 1949, p. 253). It can be shown that this approximation is justified if 


r (t) ee 
oy ta) and ae GEay PAD 


are small compared with 1. 
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The time rate of increase of the volume is therefore 


Qarr8 p(t +At) — p(t) 
SE Ax Ai (46) 


Assuming the fluid to be incompressible, this time rate of increase of the 
volume is equal to the difference between the flow into the element from 
the left-hand side per unit of time, 


awry (x), 


and the flow out of the element to the right-hand side per unit of time 


arey (x +Ax). 


Therefore we have 


arr, p (t+ At) — p(t) 


SEE x Ai =rry(x) —arr?v(x+Ax) 


or, simplifying and rearranging, 


v(%) —o(e«+Ax) 27 p(t+At) —p i) (47 
Ax FY D At : ) 


It should be realized that both the velocity v and the pressure # are func- 
tions of the two variables x and ¢. Therefore we will write (47) more com- 
pletely in the form 


v(x, t) — v(x +Ax, 2) _2r p (x, t+ At) — p(x, 2) 


Ax 6E At (ae) 
Similarly (11) should be written more completely as 
p (x,t) —p (%+ Ax, 2) v(x,t+At) — v(x, t 
oe ee ane a (49) 
Now writing (48) for ¢ + At instead of for ¢, we obtain 
v (x, t+ At) — v(a+Ax, t+ At) a 27 O(x, t+ 2At) — p(x, t+At) (50) 


Ax bE At 


Subtracting equation (48) from (50) we have, after rearrangement, 
[{ v(x, t+ At) — 9 (x, t)}—-{ 0 (x+Aa, t+ At) —0(x+Az, 2) }] /dz 


2 5 
= (FE)U ple, t+2a0 — P(e, t-+At) }—{ p(x, t+At) — p(a, t)}] /at. Og 
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Multiplying both members of (51) by p/At we find 
pl{ v(x, ¢+At) — v(x, t) }/at— {v(x-+Ax, t+ At) 


2, 
Stra Axyi) 4 ar Ae = SPY itp (x, t+ 200) (52) 


— p(x, t+At) }—{p(«,t+Ad) —p(a, ft) }]/ (an?. 


For the two terms in the numerator of the left-hand side of (52) we use 
(49) and an equation similar to (49), only with (« + Az) instead of x. 
We then obtain 


[{p (x, 1) —p (w+ Ax, ft) }/Ax—{p(x+Az, 2) 


2ré6 
OE 


—p(x,t+At) }—{p (x, t+ At) — p(x, d }1/ (As)?. 
Or, simplifying, 


—p (e+ 2s, t) }/Ax] /Ax =( [{ p(w, #+2At) (53) 


(Ax) ? 
(54) 
_27p p (a, t+ 2A) — 2p (x, t+ At) +9 (@, t) 
6E (At) 2 ; 
We see that (54) is satisfied if 
p(%+Ax, t) =p (%, t+At), (55) 
b (x +2Ax, i) =p (x, t+2At), (56) 
and 
Ax 6£ 
see ree Ba 


Equations (55) and (56) describe a wave which is propagated without de- 
formation. According to (55) the value of the pressure at the point with 
abscissa x + Ax at the time ¢ is exactly the same as that at the point with 
- abscissa x at the time ¢ + Af, and (56) has a similar meaning. The velocity 
c of propagation of the wave is given by (57) in terms of the thickness 6 
of the wall, the elasticity modulus £ of the wall, the radius 7 of the tube, 
and the density p of the fluid 
bE 

C= oat (58) 
In the case of a plus sign in (57) the wave is propagated in the direction of 
the negative x-axis, as will be obvious from the explanation above of why 
(55) and (56) describe a propagating wave. In the case of a minus sign in 
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(57) the direction of propagation of the wave is in the direction of the 
positive x-axis. Formula (58) is a special case of the formula used for the 
propagation of pressure waves in the theory of the water-hammer phenom- 
enon (Jaeger, loc. cit.), namely, for an imcompressible fluid. The most 
general form of this formula appears to have been derived for the first 
time by D. J. Korteweg (1878). Treatments taking into account the 
viscosity of the fluid (Witzig, 1914, or the Appendix to this paper) lead 
to an approximate formula in which the dependence of ¢ on 6, E, 7, and p 
is the same and only the numerical coefficient in front of the square root is 
different. 
From (58) we find 


for the velocity c in AB and 


for the velocity c’ in CD and we introduce these results into (23). We then 
obtain for the reflection coefficient y = f(t)/F(#) the formula 


9 6,4, 9 6E 
are OI: 
61Fy 6E 
R2 2 
Drip NIRS 


or, after simplification, 
(59) 


Equation (59) gives 7 in terms of the ratios 7/R, which might be called the 
coarctation ratio,” 6,/5, and E,/E. In the case of a coarctation in which 
there is only a change in diameter of the vessel 5, = 5, E, = E, and, 


therefore, 
5/2 
‘-@) 


The graph of 7 as a function of r/R is shown in Figure 4. 


For r = R, n = 1, as it should, and for r = R, n = 0, which should 
also be the case. . 


(60) 
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But if 6, # 6, and/or E, # E, or for any case in which 6,E, ~ 6E, the 
reflection coefficient 7 is given by formula (59). In this case,evenifr = R, 


n ~ 0, and the pressure waves will be reflected. 


For this case the reflection coefficient is 


aE, 

y= 24 (61) 
aE y 
6 E 


_ This would be the case if there is a local thickening of the wall, so that 
_ 6, ~ 4, as in arteriosclerosis, or other pathological conditions in which 


£, # E resulting from a locally changed structure of the wall material. 


il 


O ) 


— 
A 
R 


Ficure 4 


In all these cases the reflection coefficient might therefore give important 
clues for clinical applications. ah 
Let us now investigate how the pressure pattern might inform us about 


the reflection coefficient. Therefore we will take here as an illustration the 


_ simple case in which the incident wave is purely sinusoidal, in other words 
the functions 


a) and J (:+2) 


"are given by (25) and (26). In the latter we take ¢ = 0 which corresponds 
to m = 0 in (36). Then, as shown above in connection with (36) and (37), 


f{@ _B 
Fi) A 


4 FO 
and therefore B = 7A. 


\\ 
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The pressure at any point with abscissa x at any time ¢ is then given by 
[similar to (6)} 
p=Asinw(t—=)+n4 sino (1+). (62) 


We may write the last equation 
p= A (+n) cosw=sin wt— A (1—n)sinw —coswt. (63) 
From (63) we see that, whereas for 7 = 1 (complete reflection) 


p=2A cos w = sin wt , (64) 


which represents a standing wave with nodes for 


w= =Ftln 1=0,41,42..., (65) 
the superposition of a propagated wave and a partially reflected wave 
never produces any standing waves. In fact, the pressure represented by 
(63) has no nodes at all. This makes us wonder if in practical cases, where 
there is never complete reflection, standing waves are ever observed and 
makes us critical about claims of observed existence of standing wave 
patterns, particularly in complicated physiological cases. What can really 
be observed, at most, are minima and maxima of the amplitude of the 
wave produced by the superposition of the two waves advancing in op- 
posite directions, which occur at the so-called “pseudo” nodes and anti- 
nodes (cf. Wood, 1930, p. 281). However, in the simple case treated above, 
the latter are always a finite distance apart; in fact, a quarter of a wave 
length. However no matter what the case may be, the pattern caused by 
the superposition of an advancing wave and the wave resulting from it by 
incomplete reflection gives clues as to the reflection coefficient. As a matter 
of fact, in the simple case treated above the maximum and minimum am- 
plitudes are A(1 + 7) and A(1 — 7). By knowing A, which is determined 
by the stroke volume, 7 can be determined. From formulae (59), (60), or 
(61) we might then obtain information about 61/5, E:/E, or r/R, depend- 
ing on the circumstances. In this way clinically important new clues, 
which are of a quantitative nature, might be obtained from a study of the 
pressure pattern, if we study the reflection condition of pressure waves in 
the aorta or other blood vessels. Crude as this theory is (obviously there 
is more than enough room left for improvement) we are led to believe that 
Apéria (Joc. cit.) and Porjé (Joc. cit.) were right in stressing the develop- 


ment of a quantitative study of the reflection conditions with which they 
themselves started. 
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APPENDIX 


To derive the dependence of the velocity of propagation of pressure 
waves in a fluid enclosed within an elastic tube we will follow and extend 
the treatment given by Witzig (Joc. cit.), which is probably the most ad- 
vanced one ever published but which does not seem to have attracted 
much attention. 

The Navier-Stokes equations for the axial-symmetric flow of a viscous 
fluid in a cylindrical tube are (Lamb, 1932) 


2 
Om foes eS p (Sy 1am, Hr) 9b 66) 


Ot Wolz- a 207, Or OZ 
0 mae . ie O,\ Chas tk Os Ge ap. Op 
Ff 0 SE 3, ot) = (Sep Se Se) SP 67) 


Assuming that the gradients of the velocities are so small that the non- 
linear terms may be neglected, we have 


dv O20 ery 10 
ome 4 il a eaala 68 
Ot eerie aet+5H)- plo ox 8) 
Ov Ot e100 pee 1 Op 
—t=y¥ r r 12) pees 69 
ot DAG gE > Or me) pio ee io) 


in which y = p/p. 
Assuming furthermore that the fluid is incompressible, the continuity 
equation is: 


mee 
: 70 
z +355 an)= e feu: 
This equation is satisfied if 
1 oy 
2=—-—7z (71) 
ror 
Ld 72 
ee Pe ee 


y being Stokes’ stream function. 
Elimination of p from (68) and (69) by their differentiation to r and z, 
respectively, and afterward subtraction of the results leads to 
p eto ad wth (etees ie Seren 
la-vGatzactenca)s Gear 
From (71) and (72) we derive 
0 0, ‘ Ov, oy i Lan ay (74) 


“az. 02 rar Tor 
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Introducing (74) into (73) we obtain 


(9 (24 toy eI tse Ro 75) 


ae a2 ' ror ' or Pr r\dz2 rar 


or 
dau, (28 10 pe) (Ou eS (Aa 
jin Joteetsa a) Ae ee = r? r eke) 
For abbreviation, let us denote symbolically 
0? 1 0 0? 1 

aF oe 


027 ror art 


(77) 
Instead of (76) we then have 
(3-10) 0(8)=0. os 


Following Witzig (loc. cit.), who quotes G. G. Stokes (1819) for this ap- 
proach, we will suppose that it is possible to put 


Y=Wityr, (79) 
in which y is such that 
p(#)=0; (80) 
r 
whereas y2 satisfies 
0 Yo 
(4,-1D)@=0. (81) 
Substitution of (77) into (80) and (81) leads to 
0, Lovie Ovi 
a2 7 ayo teak (82) 


and 
A. 1 dbs | OY, 1 dye 
022 4+ Or pth way EWN 3) 


These equations allow particular solutions of the form 


Wi=F (r) e@lt-@/2)1 (84) 
and 
Yo =G (r) elt—(2/2)] , (85) 


Substitution of ¥; from (84) into (82) leads to the next equation for F(r): 


df qt dips” |e Be 
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The solution of this equation is given by 


Pir) =r jan (iSr)toni(é2 ryt, (87) 
in which 
n(itr) dl v:(i2r) 


are Bessel functions of the first order, and a and b are arbitrary constants. 
Because F'(0) is finite, whereas 


| lim Vy (i2 he 2 , (88) 
we must have 
HS ; (8 9) 
Therefore, introducing (89) into (87) we obtain 
F(r) =arJ,(i2r). (90) 
Because 


J G r) 


is real, we will write a = 7A and (90) becomes 


F(r) = Aird; (ie r). (91) 
Introducing (91) into (84) gives us 


¥i= Airs), (i2 r) eas 


(92) 
Substitution of (85) into (83) yields 
If we put hice 
Paha (94) 
we obtain from (93) Laan eee jas 
Cr ered 7, : 
which equation is similar to (86). Therefore we have 
We =BirJ, (iar) elt @/1 , (96) 
From (79), (92), and (96) we obtain 


y= r} Airs, (is r)+Birds (iar) eiatt—e/ont : (97) 


344 GEORGE KARREMAN 


To derive a relation between the stream function y and the pressure p, we 
obtain from (68), (69), (71), and (72) 


10y ya/ dp, 1dy dY\__ 14D (98) 
~ réarat roar 022 7dr OF paz’ 
1 ay 7 9 (0% _ 1,94) _19P gg) 
rosa rdz\dz22 ror dar por. 


Addition of equations (82) and (83) yields 


022 ror Or? sy Ot 


Substitution of (100) into (98) and (99) leads to 
leat ee ere 1 ap (fon 


7 drat ror y Ot paz’ 
2 
: shin 7 . a) ae -2 7 : ee 
Because of (79) we have 
and 


Multiplication of (103) and (104) by 7 and then differentiation of the 
resulting equations to z and 7, respectively, and addition of the results 
yields 

0? 10 
le de cd 2 set (105) 


02 ror ' or 


This equation allows a particular solution of the form 
p= 4H (r) elt-@/2)) , (106) 
Substitution of (106) into (105) leads to 


CH  1dH wow 
dep dp a se 


The solution of this equation is: 


H=DI.(i2r)+EY,(i2r), (108) 
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r(e2 r) snd vo(i2 r) 


are the Bessel functions of the order zero. 
Because 
lim Yo(i2 r)=o, 
r0 C 


whereas (0) and therefore H(0) are finite we must have 
E=0. 


in which 


Substitution of (110) into (108) and the result into (106) yields 
p= Je (G r) etal (¢/e)) , 


From (92), (103), and (111) we have 
D=—Apw. 


Introducing (112) into (111) we obtain 
p= —ApwJy (2 r )eialtC/0) : 
Equations (71) and (72), together with (97), yield 


v2 = 1 A = Jo (i2 r)+BaJo (iar) t eiott—(/n 


and 
-A—Ji\t4— —BaJ, iar) ¢ etl? @/1 
v7, = : Lite ‘ r aj;\1ar eA . 


Assuming that the wave length is sufficiently large so that 


wo (2r\ \* 
ic) G: S| 
NES 
ay 
p2 yp) =I, (12 are 
(i2r)= (i ee Ca 


6 if ge? 
nis rats (int) mins, 


we have, approximately, 


whereas 


and 


345 


(109) 


(110) 


(a Os 


(112) 


(123) 


(114) 


(115) 


(116) 


(117) 


346 GEORGE KARREMAN 


Consequently 
paras 4BiN (BV) t eiate— c/o, (118) 
p= — Apwel(/2)] , (119) 
v,= 1 A —+B8 Vit (1B V —i) t eiott—(/e : (120) 
n= —fid SEH BET (V9 t esol (/0 (121) 
2c¢ Cc 


These formulae were given by Witzig (Joc. cit.). 
We will now study the vibrations of the wall and take into account some 
terms that Witzig neglected. Assuming an infinitely long cylinder and 


we 


NR 
¥ 
z 
2 
oclenn nas 
~ Sie a eee 


- 
- 


= 
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axial symmetry, the equations of motion of the wall are (Love, 1944, p. 
539, Fig. 5): 


Ob: 2 Ou 

Cr Op OF ’ (122) 
aN, Td, , >, Ow 
reg ei eee oe 


in which T, and T, are the axial and azimuthal tensions, respectively, 
N, the radial shearing force in the wall, F, and F, the axial and radial 
components of externally applied forces, p the density of the wall material, 
2h the thickness of the wall, and « and w the axial and radial displace- 
ments of the considered element of the wall. 


For the axial symmetric case T,, V,, and T, are given by (Love, loc. cit., 
pp. 544, 530): 


3 ae ee 

T= DY (at oa) +p gee BH (124) 
3 2 

Ts= DI (e+ 00) —e, (125) 

N,=— one 


ae (126) 
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in which @ is Poisson’s ratio for the wall material and the flexurial rigidity 
D of the wall is given by (Love, Joc. cit., p. 503): 


2Eh3 be 
E being Young’s modulus. 
The expressions for 4, @, and k; are given by (Love, loc. cit., p. 543): 


a 
ass, (128) 
WwW 
Foe et (129) 
ow 
i (130) 


Introducing (128), (129), and (130) into (124), (125), and (126) and the 
results into (122) and (123), we obtain 


O7u Ou  o dw D2—-2¢—30@ 0w 
2ph aa = P.+5 (Satis oR M2 (ha) ore. EN 
0-w 3D /odu.w D ot2e Ow 
re SaiTAR Oia) | 2 =o) setD Ss - (132) 


The forces F, and F, are the forces exerted by the fluid in the tube on the 
wall. They are given by (Hopf, 1937, p. 95): 


Ov Ov 
ee eae E 133 
F, ee s+ or) ie 
Ov 
= hs Z : 134 
Hs p 2u(3e) ; 


after having substituted r = R in the right-hand side, as indicated. 
Introducing the expressions (133) and (134) into (131) and (132) we 
have: 


02u oer 0 0; o Ow 
2 pha = — # ats +5 a Rdz 
(135) 
D 2—20—30? dw 
RO? (Lh 6) 023 | 
ow 0 Ur _3D/(a ou 
2ph— =p—2p a ot 


D ot2o? dw 7 dt 
~R2(1—«a) 02? d24° 
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Assuming no slip of the fluid at the wall, we have also: 


Ou 


er Grae (137) 
Ow 
ee (138) 


Differentiation of equations (135) and (136) to ¢ and introduction of the 
results (137) and (138) gives us, after neglecting the terms which contain 
03 0, 0490, 


Veet er 


(which are supposedly sufficiently small), for r = R: 


Oot 0, 07 0, Ov, 6 00; 
OP genes sate) az? +35) (om 
C25 ene p) 0" 0, Ge . Oh, © GEG 
bf? ais aiar NR Oe oe (140 
with 
preg 
34D) E 
2 ph® pA — 0)’ Le 
Sas 143 
ome ye ( ) 
ED see 2s" 
Sr naa (144) 


For y = 0 and ¢ = 0 (139) and (140) reduce to Lamb’s equations [cf. 
Witzig (loc. cit.)). 


From (119), (120), and (121) we have, for r = R, 


= — A pwelt(2/e)] , 


= 44 “+BB ViJy (RB V1) telat ceron, (145) 
: 7 
= — {id oe +B 2 71 (RB Vi) f ciate) (146) 


Introducing these equations, together with (145) and (146), into (139) 
and (140), leads to 


aA+yB=0, (147) 
KAS AB = 0; (148) 
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with 
w? 1 1 
a= )1—S aR 5-14 (1-3) . (149) 
~ Jo (RB -V —i) 7 
= oi RBV— —j Na (1 
iF Ji (BBV =i) {RB V— ia Eh (1 5) os 
MESS apse e. - 
(6) c* Cc 
c= sR ibpw ty S +5 $2412" (4), (151) 
n= 2, (RBV =i) [ a? 2 Hh _ tis Voi pea) 
£ = Jo (RB V — =i) ae 
ww” mae A ari 
BR pe kbV— CREA a) ak 
Elimination of A and B from (147) and (148) leads to 
at=yYk, (153) 


or, after introducing (149), (150), (151), and (152) into (153) and sim- 
plying, 


ee Shi ie ? Seth 1+RBV—i 
swS) pe. ee rae ae | a 
Bon ae ae 
to BES Rien ty S45 pati Fa (ED |. 


Neglecting the inertia terms and some other terms which are negligible for 
plausible values of the parameters: y~1, R~1, w~ 10, c ~ 10%, 
7 ~ 10°, ¢ ~ 103, o ~ 0.5, and dropping some terms which Witzig (Joc. 
cit.) did not take into account, we find as a check Witzig’s formula 
THEME Ly PO) 155 
{Re (c) iP Rp Rp(P—o)?+C ’ ( ) 
in which Re(c) is the real part of c, and P and Q are the real and imaginary 
parts, respectively, of the function 
- Jo (RB -V =i) 
RBV i + (RBA 4) 


which occurs in (154). 
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Still neglecting the inertia terms, but keeping the terms 
2 , 4 
TRE and) ee 
@ 2 Ge 
which are lacking in Witzig’s treatment but which may not be negligible, 


7 taj 
we obtain _Eh (P — 2) (P—4) +e 


Rp @-F?+@ 
For ¢ = 4 we again obtain Witzig’s formula. Keeping both inertia terms 
and the terms mentioned above, we obtain after some lengthy numerical 
calculations the next result: 


Re(c) = VEE (0.8 + 0.03Ry2). (157) 


As expected, the influence of the viscosity on the velocity of propagation is 
very slight. The damping is given by the dependence of Im/(c) on y which 
can be calculated in the same way. 

The author would like to acknowledge the assistance of Mr. D. L. Cohn 
who made the calculations which led to formula (157). 

This work was made possible by a grant from the National Heart Insti- 
tute to the University of Chicago. 


{Re (c) }? (156) 
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The response time of a random net is defined as the expected time (measured in 
the number of synaptic delays) required for the excitation in the net (measured by 
the fraction of neurons firing per unit time) to reach a certain level. The response 
time is calculated in terms of the net parameters as a function of the intensity of the 
outside stimulation. Two principal types of cases are studied, 1) an instantaneous 
initial stimulation, and 2) continuously applied stimulation. It is shown that for a cer- 
tain type of net where the required level of excitation is small, the response time- 
intensity equation reduces to the one derived on the basis of the “‘one-factor’’ theory 
applied to a neural connection. More general assumptions, however, give different 
types of equations. 

The concept of the ‘‘net threshold’’ is defined, and its calculation indicated. The 
net threshold for instantaneous stimulation is, in general, greater than that for con- 
tinuous stimulation. The results are discussed with reference to existing theories of 
reaction times. 


The treatment of response times as a function of stimulation intensity 
appears in mathematical biology in two connections; 1) the response of a 
peripheral nerve to an applied potential difference, which may be constant 
(Blair, 1932a) or variable, as for example in a condenser discharge or in a 
source of an alternating current (Blair, 1932b); and 2) the response of the 
total organism to stimulation of its sense receptors, as in the determination 
of reaction times (Piéron, 1920-1921; Householder and Landahl, 1945; 
Rashevsky, 1948). 

In the case of the peripheral nerve the simplest mathematical treatment 
of response time-intensity relations is that of H. A. Blair, sometimes re- 
ferred to as the ‘‘one-factor theory.” There it is assumed that the “‘ex- 
citatory factor,” whatever its nature might be, is subjected to two influ- 
ences, namely a tendency to accumulate at a rate proportional to the mag- 
nitude of the impressed voltage (or the current flowing through the nerve, 
if the resistance of the latter is assumed constant) and a tendency to dis- 
sipate at a rate proportional to its own “concentration.”’ Hence the 
equation 
—=Al—ke, (1) 
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where ¢ is the “concentration” of the excitation factor, 7, the impressed 
voltage (or current), A and k being constants. : 

It is further assumed that the nerve becomes active, i.e., the spike po- 
tential appears and begins to propagate along the nerve, when the concen- 
tration of excitation reaches a certain critical level /. Integration of equa- 
tion (1) and the substitution of / for e then gives the response time ¢* as 
a function of the intensity J and the constants (the parameters of the 


nerve), namely 
) 1 Al 


2 Saar 


je (2) 
which satisfies the conditions é*(h) = , f*() = 0. Equation (2) gives 
good agreement with experimental data on the response times to a con- 
stant potential applied to peripheral nerve (Blair, 1932a). 

N. Rashevsky (1948) has used a formally equivalent theory to derive 
the response time of an efferent neuroelement as a function of the intensity 
of activity in an afferent neuroelement synapsing on the former. It was 
found, furthermore, that if the intensity of activity in the afferent neuro- 
element was assumed to be a logarithmic function of the external stimulus 
intensity (e.g., brightness in the case of visual stimulation, concentration 
in the case of a gustatory stimulus, etc.) very good agreement could be ob- 
tained between the corresponding theoretical equation, relating stimulus 
intensity to response time, and the reaction times of human subjects to 
sense stimuli (Householder and Landahl, loc. cit.; Rashevsky, loc. cit.). 

It thus appears that the response time-intensity relation derived from 
the simple assumptions of the one-factor theory (plus a logarithmic inter- 
pretation of stimulus intensity) fits the data on reaction times, where 
probably large sections of the nervous system are involved. It seems at 
first surprising that this should be so, inasmuch as only one neural connec- 
tion is involved in the model where the equation of the one-factor theory 
is applied to the derivation of the over-all reaction time; all the more so 
since the magnitudes of these reaction times exceed by a factor of some 
hundreds the usually observed magnitudes of synaptic delays. N. Rashev- 
sky (1948) has shown, however, in his reinterpretation of the fundamental 
equation in terms of large numbers of reverberating cycles, that the “single 
connection” can be interpreted as a “macro-synapse” with a large number 
of neurons, so organized that the continuous rise of excitation is a reflec- 
tion of the activation of more and more cycles, while the dissipation of ex- 
citation is due to accidental de-activation of the cycles. The formal aspect 
of the theory is thereby preserved. 


In this paper we shall propose an alternative model, also based on a 
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large number of neurons, each behaving in accordance with the all-or-none 
law, where no assumptions concerning the organization of these neurons 
into cycles are required. On the contrary, we shall suppose the greatest 
possible chaos, characteristic of a random net, where the probability of a 
direct connection between any pair of neurons is the same. Such a random 
net was described in previous papers (Solomonoff and Rapoport, 1951; 
Rapoport, 1951, 1952). We shall refer most often to the last mentioned 
_ paper, hereafter designated by Joc. cit. 

If, as it is reasonable to suppose, the neural events involved in reaction 
times of human subjects following instructions (i.e., events other than 
simple reflexes) depend on the activity of higher centers (e.g., the associ- 
ation areas of the cortex), it seems to the author advantageous to use prob- 
abilistic concepts in developing a theory of such events. For, loosely speak- 
ing, it can be conjectured that such centers are much less “rigidly” organ- 
ized than the older pathways involving, say, inherent reflex responses. The 
functional pathways of the higher centers may well be in a state of relative 
flux, quickly changing with the situation to be met. That is why the pos- 
tulate of a chaotic random net seems to the author a more promising ap- 
proach than postulates of specific connections, specifically organized to 
explain a psychological phenomenon. 

Accordingly we shall consider the response time of a random net of the 
type described in previous papers, particularly in loc. cit. Some of the 
neurons of such a net Yt will be supposed to be stimulated from the out- 
side, either instantaneously (just once) or continuously (once every unit 
of time). Quantized time is assumed throughout, the unit being the syn- 
aptic delay. Any imposed excitation may spread throughout it because of 
the connections between the neurons of St. The response time will be de- 
fined as the time required for the activity of Jt to reach a certain level, as- 
sumed as being the threshold of the effectors stimulated by 9t, from which 
effectors the behavior associated with the overt reaction to the stimulus 
results. 

In loc. cit. we designated by x(t) the number of neurons firing in Jt at 
the (quantized) instant ¢. In this paper we shall consider instead a(t) = 
‘2/N, i.e., the fraction of neurons firing. The translation of the equations of 
loc. cit. into this notation is obvious. 

Case 1. Instantaneous Stimulation; h = 1; a > J. It will be recalled (cf. 
loc. cit.) that h is the individual threshold of the neurons of %, ie., the 
number of stimuli which must impinge simultaneously upon any neuron 
in % to elicit a firing. If each stimulus received by a neuron results in a 
firing, then 4 = 1. The axon density, denoted by a, is the average number 
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of axons sent out per neuron. In the case we are about to consider, we sup- 
pose that k = 1, a > 1 (the case # = 1, a S 1 is trivial, since any initial 
stimulation will die out, i.e., z = 0 is the only steady state). As stated 
above, z(é) is the fraction of neurons firing at the instant ¢. Following the 
reasoning of Joc. cit., we arrive at the equation 

dz 

—=1-—e%—2z. (3) 

di : 

If H is the level of excitation at which response takes place™ [i.e., the 
“over-all” threshold of the net, H > 2(0)], then the time of response will 
be given by 

H 
t* = BASES aol (4) 
Hoe na Opetane 
Equation (4) gives the response time-intensity curve under the conditions 
considered. Here z(0), being the initial and only outside stimulation, rep- 


x ©) 


Ficure 1. Response-time plotted against instantaneous stimulation for various values 
of a. The lowermost curve is the limiting case, a = ~. The dotted line is the undefined case, 
a=1. 


resents the intensity, while H and a are the parameters of the net. By in- 
spection, we see that t*(0) = ~; ¢*(H) = 0; dt*/dz(0) < 0 for all values 


0 <2(0) < H, (a >1). For very large a (negligible e-), we have the 
approximation 


=loe =, AG 


while for a very close to 1, equation (4) becomes meaningless since the 
steady state approaches zero and H cannot exceed it. The family of curves 
for our range of a is represented in Figure 1. 


* We suppose that H is less than the steady state level of activity as defined in loc. cit. 
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If aH <1, it is easy to obtain an approximate quadrature indicated in 
(4) in closed form. In that case az < 1, and we can approximate e- by 
1 — az, obtaining 

dz 1 


H 
aw = et 
eo 21) woe a0) - se 
With a better approximation, 
2 o2 
é V@ail—azt “S ; 
formula (6) is corrected to 
= ne 
Oo n(@atpe 
t* ~~——__ log ; "ae (7) 
a—1 = 


Any desired approximation can be obtained by expanding the integrand 
of (4) in series. We first expand the denominator getting 


1— e-*%—z=2(a—1) apes (1) 7 (az)? 
7=2 (8) 


| (a— 1) 4 Soe 


j=l 

The series within the brackets of (8), being analytic at z = 0, can be in- 
verted. Its reciprocal is 

1 a a!+2a3 


2+ 


PRISON = Tee TT (9) 


Hence the entire reciprocal of the integrand is 


1 a at+ 2a 


he eA 0 
Ga? ite nis ’ (10) 


which, integrated termwise, gives us the expression for ¢* as a series in 
2(0) and H, namely 
Hf a 


1 Ae ee We RE 
a sat 8H Tage zi 
Se LHP = 22 (0)1 + : set 


This series converges quite rapidly if a is not too large, i.e., aH < 1. 
Case 2. Continuous Stimulation. Now let f be the fraction of neurons per 

unit time in 9t constantly stimulated by an outside source. Hence f is a 

measure of the intensity of that stimulation. Pick a neuron of 9 at ran- 
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dom. Consider the probability that at time ¢ + 1 it is stimulated neither 
by the outside source nor by any of the neurons within 9 [of which there 
were x(t) active at time f, ie., having ax(¢) = Naz(é) active axones]. This 
probability is given by the product 


(1-7 \ccmne (11) 


The first factor is the probability of not being stimulated by the outside 
source; the second, the probability of not being stimulated by a neuron of 
N. The two events are assumed independent. Therefore the probability of 
being stimulated by either source at the time ¢ + 1 will be 


De lea) Crees (12) 
and our differential equation for the spread of excitation becomes 
—-=1— (1—f)e*—s; (13) 


The reaction time will accordingly be given by 


H dz 
t= f re ee (14) 


Of particular interest is the case when a < 1, H <1. We recall that this 
case was trivial under the condition of instantaneous stimulation, because 
there the excitation simply died out, whatever was the outside stimula- 
tion. Here, however, excitation is constantly “fed” into our net. If H (and 
consequently the range of z) is so small that all but the linear terms may be 
neglected, equation (14) becomes upon integration 


= 7a log (15) 
which is of exactly the same form as the “classical” equation (2) with 
AI = f (the activity of the afferent); 1 — a = & (the rate of ‘“‘dissipation” 
of the excitation) and H = h (the threshold of the efferent). 

With a proper interpretation of f as a function of the intensity of the 
outside stimulus, we have here a random net model of the reaction-time 
theory. 

The condition a < 1 becomes particularly interesting if one allows one- 
self a little speculation. In our random net, a is the axon density of the net. 
Small a (a < 1) means that on the average less than one axon is emitted 
per neuron. Anatomically speaking, this assumption is quite improbable, 
since it postulates the existence of neurons without axons. Functionally 
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speaking, however, the interpretation is different. Our “net” is not neces- 
sarily a geographically determined collection of neurons. It is rather 2 
functionally determined collection which becomes functionally significant 
in a given situation (for example the subject’s attention being fixed on the 
requirements of the experiment). Hence our net may be a selected collec- 
tion of neurons which will participate in a given task. In such a collection, 
the axon density may well be less than unity. With practice, certain neu- 
rons may be functionally eliminated from the collection and certain others 
added with a resulting increase in the axon density and a corresponding 
reduction in reaction time. We note from Figure 1 that this reduction is 
most pronounced (in the case of instantaneous stimulation) at small 
values of z(0), indicating greater “sensitivity” to proper stimulation. Simi- 
lar effects may be derived from the assumption of continuous stimulation. 

Case 3.h > 1. We have seen (Joc. cit.) that, for k > 1 and under the as- 
sumption of quantized time, excitation in a random net started by an in- 
stantaneous stimulation propagates according to the following equation: 


d 
sat- e “#,_, (az) —2 (16) 


where 
k 
yi 
EF; (y) = 2 re . 
rarten 


The corresponding response time will be 


az 


ope ead 
Sg) bre hes (a2) 


(17) 


We have also seen that random nets with # > 1 exhibit the “ignition 
phenomenon,” that is, the initial excitation of a certain minimum number 
of neurons will result in a further spread of the excitation toward a stable 
steady state, whereas the excitation of a number of neurons less than the 
minimum will die out. This “ignition phenomenon” is reminiscent of a 
threshold. However, the notion of threshold in peripheral nerve excitation 
is usually connected with a minimum applied voltage sufficient to excite. 
Analogous to this notion would be the continuous stimulation of a con- 
stant fraction of neurons. ' 

Let f be that fraction. Then, following the same reasoning as that pre- 
ceding equation (9), we obtain the modified form of equation (16), namely 


da _ P23 = ¢- Yi (02) <2. (18) 


dt 
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The proof of the existence of the ignition phenomenon in loc. cit. pro- 
ceded by a semi-quantitative argument. It will be recalled that we super- 
imposed the graphs of 1 — z and of 


eH ,-1 (a2) 
against z over each other, and showed that, for sufficiently large values of 
a, the line 1 — z intersects the curve 

e~E,—1 (a 2) 


in two points, z = 6: and z = 72, the two steady states of z. We further 
showed that 71, the smaller value, was an unstable state, and y2 a stable 


(pe Oe?) 3 


oR 1-3 


Prep (a3) 


% z, 4 
Ficure 2. The line 1 — z and the family 
(1 — f)e-* Ey (a2) 


plotted against z for various values of f. The uppermost curve is for f = 0. The next curve is 
for f < H. The tangent curve is for f = H, and lowermost for f > H. 


one, so that the excitation starting with z < y, would tend to die out, 
while one starting with z > y, would tend toward 72. The reader is referred 
to Figure 1 of loc. cit. 

Let us now apply similar considerations to the right side of (18). We 
hold a fixed at a sufficiently large value so that for f = 0 there are two in- 
tersections. Consider now a family of curves corresponding to increasing 
values of the parameter f, which, we recall, is a measure of the intensity 
of the outside stimulus. This family is shown in Figure 2. 
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As we increase f, the factor (1 — f) of the expression 
Cl = jrer he (a2) 


decreases. Thus the entire curve is depressed. There is a critical value 
f* = H for which the curve is tangent to the line 1 — zg. This critical value, 
we maintain, is the threshold of the net % under continuous stimulation. 

To establish this, we must show that for f < H, the initial fraction of 
neurons stimulated is less than +; (the unstable steady state) and remains 
so, otherwise the excitation will tend toward y2. That this is so follows 
from the fact that the derivative with respect to z of 


Cea 1 Gas) 


is always negative for z > 0, as can be verified by direct differentiation. 

This implies that the ordinate of any point on the curve is less than the 
ordinate of the starting point (where z = 0). If f < H, there are three in- 
tersections of the curve with the line 1 — z corresponding to values Yo, 1, 
and ‘2 of z. Therefore the ordinate of yo is 1 — yo, and this is less than the 
ordinate at z = 0, namely 1 — f. Then 


1-—-y<1i-f or i KN - (19) 


Hence, a forteriori, f < 1, and our contention is proven with respect to the 
initial stimulation. That excitation in this case will tend toward yo < 71 
follows from the fact that at z < Yo, dz/dt > O and at z > yo, dz/dt < 0. 
Therefore yo is a stable state. Therefore, for any subthreshold stimulation 
f < H, excitation will at first grow until it is stabilized about yo, a nonzero 
steady state. 

For f = H, yo = v1 becomes an unstable steady state, and the excita- 
tion will spread to the next stable steady state y2. For any value f > H, 
there is only one steady state, namely, 2. The evaluation of the threshold 
H can be obtained from the following 

Theorem. Leth > 1 be the individual neuron threshold in a random net Nt, 
whose axon density is sufficiently large to give non-zero steady states. Then 
the over-all net threshold H with respect to a constant stimulation f 1s given by 
the simultaneous solution of the following equations: 


f= (1—/f) ¢ °#,-1(¢2) —2=90, (20) 


qh-l gh-l 


a(1—f) Sr esa Wa ad (2:1) 


Proof. Equation (20) determines the intersections of the curve 


(1 — f) .e—*#Hy,-1 (a2) 
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with the line 1 — z, and hence the steady states. Equation (21) is obtained 
by setting the derivative of 


(1 — f) e~*B,-1 (22) 


equal to —1. The two equations together are the conditions for tangency 
between the curve and the line, in which case yo = 1, ie., f = H. 

Remark. It must be kept in mind that all our considerations are proba- 
bilistic. It is not certain that a subthreshold stimulation results in the ex- 
citation always staying below yy. This value is only a theoretically de- 
duced steady state. Actual excitation may fluctuate around it, being some- 
times above and sometimes below. Now if yo and 1 are close together, i.e., 
fis subthreshold but close to the threshold, it may happen that a fluctua- 
tion of g may carry it beyond 1, and an “explosion” toward 72 may result. 
We have here a model for over-all threshold fluctuations, even though the 
individual neuron thresholds are fixed at k = 1. It would be interesting to 
compute the form of these fluctuations, for example, to determine if they 
are normally distributed. 

Conclusions. We have shown how a theory of reaction times can be de- 
rived on the basis of a “discontinuous” theory of neural action without 
postulating anything about the structure of an aggregate of neurons other 
than the condition of ‘‘chaos.”’ 

In a very special case, the equation deduced from the random net theory 
reduces to the “classical” equation derived from the application of the 
one-factor theory to a single connection. Hence the random net equations 
are a generalization of the one-factor theory equations. However, the most 
interesting theoretical result of the random net theory of reaction times is 
that it implies a quantitative relation between the thresholds of stimula- 
tion for two methods of stimulation—by a short (flash) stimulus and by a 
continuous (constant) stimulus. The former threshold can be computed 
from equations (20) and (21). The latter threshold is the smaller root of 
equation (21) of loc. cit. For fixed parameters of a given net (i.e., a and h), 
it is easily seen that the threshold of continuous stimulation is always 
smaller than the threshold of flash stimulation. This result finds qualita- 
tive experimental verification (Piéron, loc. cit., Geblewitz, 1935). 

Furthermore a quantitative relation is implied between the respective 
reaction times for the cases of super threshold flash stimulations and con- 
tinuous stimulations of equal intensity. The former is given by equation 
(17); the latter can be obtained analogously from equation (18). Even if 
the continuous stimulation is a given function of time f(t), the reaction 
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time can, in principle, be calculated from equation ( 18). These additional 
relations provide a possibility of reducing the number of unknown parame- 
ters in the reaction time equation. 

One is tempted to apply these theoretical results to available data on 
reaction times. A wealth of such data, especially relating to the compari- 
son of reaction times for continuous and flash stimulation is found in the 
work of Piéron and his school (Piéron, Joc. cit.). It appears, however, that 
such applications cannot be made without a rather thoroughgoing 
analysis. 

To begin with, Piéron’s scale of stimulation intensities is not in terms of 
some objective unit but in multiples of threshold stimulation. Since the 
thresholds for the two methods of stimulation are widely different, the 
data must be reinterpreted in objective (energy) units. 

The second difficulty is more serious. The wide variations of reaction 
times for various modalities of sensation indicate that perhaps very differ- 
ent classes of events are operating in each case. What is measured in any 
reaction time experiment is the time interval between the activation of the 
stimulation source and a previously agreed upon response in the subject 
(in the case of human subjects, the pressing of a key). During this time, 
the following events presumably occur: 1) the passage of energy from the 
source to the receptors of the subject; 2) the activation of the receptors; 
3) the passage of impulses along the neural pathways of the subject; 
4) the contraction of appropriate muscles which constitute the response. 

Strictly speaking, a theory which deduces a quantitative relation be- 
tween the intensity of the stimulus and reaction time must relate the in- 
tensity to the time course of events in each of the above mentioned phases. 
The usual practice is to lump three of those levels into a constant and to 
deduce the intensity-reaction time relation in terms of the remaining one. 
This is implied in one of the earliest theories of reaction time, namely, that 
deducible from the so-called Bloch’s Law, which postulates a simple re- 
ciprocal relation between intensity and the time necessary to excite (con- 
stant “dose’’). Here the variable factors are by implication attributed to 
the events at the receptors, all the other events (including those in the cen- 
tral nervous system) being considered constant. 

On the other hand the Landahl theory (Householder and Landahl, loc. 
cit.) considers the variable component to be the events at a connection be- 
tween an afferent and an efferent neural pathway and lumps into a con- 
stant the time passing in the activation of the receptors and in what occurs 
after the excitation has crossed the connection. 

Piéron (Joc. cit.) attempts to consider both the activation of the recep- 
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tors and the neural events as variable, but the latter are not treated with 
mathematical rigor for want of a suitable theoretical model. 

It appears that three rather different theories yield similarly shaped 
curves relating intensity of stimulation to reaction time. Piéron’s theory 
is essentially a modification of Bloch’s Law. It focuses attention on the 
activation time of the receptors. The Landahl theory is an application of 
Blair’s one-factor theory to a neural connection together with an assump- 
tion of a logarithmic relation between the objective intensity of the stimu- 
lus and the activity in a neural pathway. The random net theory assumes 
a spread of excitation through an aggregate of neurons randomly connect- 
ed and interacting “digitally,” i.e., in accordance with the discontinuous 
model of the central nervous system of McCulloch and Pitts (1943). The 
similarity of the intensity-reaction time curves deduced from the three 
theories makes it possible to fit the observed data to any one of them. To 
compare the relative merits of the theories, one must therefore go beyond 
curve fitting—one must find suitable interpretations of the parameters in- 
volved and devise means of varying these parameters in an experimental 
situation. This leads to questions concerning the relative importance of the 
different phases of activity listed above in the interpretation of the reac- 
tion time curve. It is noteworthy, for example, that the reaction times for 
certain modalities of stimulation are much larger than others. For ex- 
ample, in thermal stimulation of the human skin (Geblewitz, loc. cit.) the 
response for continuous stimulation are of the order of one second, ap- 
proaching asymptotically the minimal value of 0.3 to 0.4 seconds. For 
gustatory stimulation (Piéron, Joc. cit.) the response times in the vicinity 
of the threshold are of the order of 2 to 3 seconds, approaching the mini- 
mal value of 0.25 seconds. Experiments on invertebrates have shown re- 
sponses to light stimuli delayed as long as 7 seconds (Hecht, loc. cit.). Re- 
sponses to visual stimuli, on the other hand, are usually in the vicinity of 
0.3 seconds in the threshold neighborhood and approach values of less 
than 0.2 minimally (Piéron, Joc. cit.). 

These wide variations indicate that many different mechanisms may be 
responsible for the intensity-response time relation. Any refinement of the 
existing theories, therefore, must take some such possible mechanisms 
into account. In the case of thermal stimulation, it would be interesting 
to consider the quantitative relation between the intensity of stimulation 
and the time necessary for the energy of stimulation to penetrate the epi- 
dermis to the receptors. In the case of gustatory stimulation, the same 
problem presents itself in connection with the diffusion of the chemical 
stimulant throughout the region where the receptors are situated. This 
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latter problem is especially interesting inasmuch as different chemicals 
give different reaction time curves. 


Some of the considerations outlined here will be treated in subsequent 
papers. 


This work was aided by a grant from the Dr. Wallace C. and Clara A. 
Abbott Memorial Fund of the University of Chicago. 


LITERATURE 


Blair, H. A. 1932a. ‘“‘On the Intensity-Time Relations for the Stimulation by Electric Cur- 
rents, I.’’ Jour. Gen. Physiol., 15, 706-29. 

. 1932b. “On the Intensity-Time Relations for Stimulation by Electric Currents, II.’’ 
Ibid., 15, 731-55. 

Geblewitz, E. 1935. ‘‘La relation entre l’intensité du stimulus et le temps de réaction pour 
les excitations thermiques en durée indéfinie et en durée bréve.’”’ C. R. Soc. de Biologie, 
119, 1362-1364. 

Hecht, S. 1918. ‘‘The Photic Sensitivity of Ciona Intestinalis.’’ Jour. Gen Physiol., 1, 147-166. 

Householder, A. S. and H. D. Landahl. 1945. Mathematical Biophysics of the Central Nervous 
System. Bloomington, Ind.: Principia Press. 

Piéron, H. 1920-1921. ‘“Nouvelles recherches sur l’analyse du temps de latence sensorielle et 
sur la loi qui relie ce temps 4 l’intensité de l’excitation.’’ L’ Année Psychologique, 22, 58-142. 

Rapoport, A. 1951. ‘“‘Nets with Distance Bias.’’ Bull. Math. Biophysics, 13, 85-91. 

. 1952. ‘‘Ignition Phenomena in Random Nets.’ Ibid., 14, 35-44. 

Rashevsky, N. 1948. Mathematical Biophysics. Rev. Ed. Chicago: University of Chicago Press. 

Solomonoff, R. and A. Rapoport. 1951. ‘“Connectivity of Random Nets.”? Bull. Math. Bio- 
physics, 13, 107-17. 


BULLETIN OF 
MATHEMATICAL BIOPHYSICS 
VOLUME 14, 1952 


THE SMALLEST VALUE OF THE AXON DENSITY FOR 
WHICH ‘IGNITION’ CAN OCCUR IN 
A RANDOM NET 


ERNESTO TRUCCO 
COMMITTEE ON MATHEMATICAL BIOLOGY 
THE UNIVERSITY OF CHICAGO 


As shown by A. Rapoport (1952), when a very brief stimulation or ‘‘nstantane- 
ous input’’ is applied to a random net, the subsequent events are determined by the 
parameters of the net as follows: If the axon density a is sufficiently large and the 
fraction y of the neurons initially stimulated exceeds a certain value y; (the over-all 
threshold of the net for instantaneous stimulation), excitation will spread through 
the net until a steady state is reached in which a fraction y2 > v1 of the neurons is 
firing (‘ignition phenomenon’’). If y < +; the activity in the net dies out. However, 
if the axon density is too small, the activity will ultimately die out, no matter how 
large the fraction of initially stimulated neurons. Thus there exists a limiting value 
A of the axon density below which the net cannot “‘ignite.’’ This A is a function of h, 
the individual threshold of the neurons constituting the net (we assume here h > 2, 
since for k=1 the situation is essentially different). Geometrically y1 and ye are de- 
termined as the two intersection points of a straight line with a sigmoid curve. When 
a <A the two curves do not intersect and for a = A they are tangent. 

In this paper the ‘‘tangency case’’ is investigated and the general features of 
the function A(z) are determined. It is shown that A increases monotonically with 
h (as one would expect). For all values of > 1 we have A(h) > h, but the fraction 
A(h)/h and the derivative dA(hk)/dk approach unity as # increases. An analytical 
expression of the function A(/) valid for very large values of & is derived. 


Consider a random net Jt as described by Rapoport (oc. cit.). Upon be- 
ing fed an initial instantaneous input, the net can reach a positive steady 
state only if 

cach ay) y= 0) (1) 


[Rapoport, equation (20)]. 

Here + is the fraction of neurons firing steadily, a is the axon density 
(average number of axons emitted by each neuron) and # (an integer) the 
individual threshold, assumed equal for all neurons. 

Fora given k > 1, equation (1) has no real solution y if ais smaller than 
a certain number A (which we wish to determine), whereas if a > A equa- 
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tion (1) has two distinct solutions yi and 72. This is seen by plotting the 


two functions 
OAY)) elegy, 


and 


ll 


v (y) e~ °1H,-1 (a7) 


vs. y (Fig. 1; also Fig. 1 in Rapoport). 
For a = A the two curves ¢(y) and w(y) are tangent, i.e., we have: 


1—-y=e %E,-1(AyY), (2) 
a AQ 
Sle —A¢ Tea i: (3) 


By eliminating y from (2) and (3) we obtain a relation between A and 


FIcuRE 1 


h, giving us the minimum value of the axon density for which the net will 
still “ignite” as a function of the individual] threshold. 
If we put 
Aq y (4) 


equations (2) and (3) become: 
l1—-y> =e F,-1 (9); (S) 


—=€e_ 


Shere err (6) 


Inserting (6) into (5) we obtain an equation involving only f/ and y: 
= ey 
Le ay) mie cae (7) 


This last equation implicitly defines a function y = f(z), which, when sub- 
stituted into (6), will show how A depends on h. 

It is easy to solve equation (7) by trial and error for small values of h. 
We have, e.g., for h = 2: 


Lei ey reves 
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which is satisfied for y = 1.79... = y(2). We then find from (6) that 
Pl) = 3.34... 


Similarly, for k = 3 equation (7) becomes 
i e+ 1 +y+2| eee 
Y ? 


which gives us y(3) = 3.38... and A(3) = 5.14.... 
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Ficure 2. Ordinates at right refer to dash-lined curves, ordinates at left to full-drawn 
curves. 


Figure 2 shows the graphs of y(#) and A (A) calculated in this way for 
values of # ranging between 2 and 11. (Numerical data are collected in 
Table I.) 

However, to obtain a complete picture of the function A (/) it is neces- 
sary to study the behavior of equation (7) for very large values of h. 
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Putting 
G,(y) =1— eH (y),; (8) 


we notice that 


1 y 
G;, (y) Tea | if e tt —ldt . (9) 


Equation (9) can be shown by induction or, more directly, by verifying 
that 


qGa(y) _ op Soo 
dv (h—1)! 


and G;, (0) == (0). 


TABLE I 
5 A(h) calcu- 
h y(h=n A(h) n—h | (t—1)/n | 1/A=y_ | lated from 
equation (50) 
PR aS 1.79 3.34 | —0.21 0.56 O54: ail cep onions 
Simnaee 3.38 5.14 0.38 0.59 (UR esc Oe eset mr: 
As aes 4.88 6.80 0.88 0.61 OST 24s ae ae 
avd hero 6.32 8.36 Wao ys 0.63 OG 2" Pea 
One aoe thst? 9.86 1 2s 0.65 0578. 394-2 ee 
1 tees 10.45 12.79 2.45 0.67 0.82 9.47 
1 eee 14.39 16.93 3.39 0.69 0.85 13.67 
Equation (7) now becomes: 
ev y Pets 
mee Midt=1, (10) 


or, introducing the new variable of integration, u = 1 — (¢/y), 


Heo =f enw det. (11) 


It may be of interest for future calculations to notice that the integral 
in (10) is the incomplete T'-function y(h, y) (in the notation of Whittaker 
& Watson, 1945, p. 341), whereas for the integral in (11) we have: 


Wy) =F-F UL b+) 591, | (12) 
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where 


mit, th 1); y| Sih ae rea 


2 


+1) (h+2) 
+ 


We 
i GG ee 


is the well-known confluent hypergeometric series.* 
Equation (7) has only one solution y(k) ¥ 0 (y = 0 being a trivial solu- 
tion). This can most easily be seen from (11). We have: 


dH, (y) 


1 
iy =f ev (1—u)*"Iudu>0. (13) 
0 


Thus, H;(y) and dH,(y)/dy are monotonically increasing functions of y. 
Also: 


H, (0) =7<! Te Rae) 


dH, (y) 1 oe 
ny ne 
| dy ae ae mate 


and, for any fixed value of y, 
Anis(y) <Any), } 


dHnii(y) — dHn Cy) f (15) 
dy Srey, 


Plotting the family of curves H;(y) vs. y (Fig. 3), we therefore see that 
any of them intersects the horizontal line of height 1 in one point only. 
Moreover, (2) increases monotonically with h. 

We also plot the functions G;(y), (8), and 


nate 
M;, (v) = eaA)t ! (16) 
vs. y (Fig. 4). Equation (10) can be written: 
Gr(y) = Mi Cy) - (17) 


* In general: 


PCY) LNG) ay 
F(a, 73 9) Soran vt 


= TG Fr 1 ev (1—u)121u“"1 du, 


F(a, 7; 9) =1, 


LaF (a, 73 9) i = 


a 
(iy Ge 
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Now, G,(y) remains less than 1 and tends to 1 for y—> , whereas 
M_,(y) has a maximum at y = h and becomes equal to 0 for large y’s. 
Using Stirling’s approximation, 

ia 
hyhow (h— 1)! Ve— 18 
e hhh Ne, (18) 


we find: 


M, (h) ae. (19) 


FIGURE 3 


y=h yih=7 y 


Ficure 4. Not on scale 


This shows that 


; y (h) >h (20) 
if h is sufficiently large.* 


In what follows we shall denote by 7 [=y(h)] that particular value of 
y which satisfies equation (7) for a given h. Then, multiplying (7) by 
exp(n — /) and using (19) and (20), we obtain: 
a _ fe oy oe ee 
ait ta 
* Actually (20) holds for h > 2 (see Fig. 2). 
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Hence: 
n—h>¢$ln a), (BAN 
lr Y 
n 1 h 


Since (i) increases with 4 and G,(y) tends to 1 for y > ~, we would 
expect that for large h’s equation (10) or (17) should be equivalent to 
ae eeaBE nS 
OS ampa Qs) 
This can be shown rigorously as follows. Multiplying equation (7) by 
[1/M;(n)] we find: 


ees) 
é€ ™7 
where 
— { 
Cee a Sener. 


For an arbitrarily large but fixed # we then have, since (hk — 1) /n < 1, 


a ypu ah 1) ee 
b(n) = S Se 


ane n 
qh et lg et (h—1— yp)! 


te al Ne 1 
LC eee 
er aera 


Using (21) we therefore conclude that ®(n) becomes negligibly small with 
increasing / and this proves our assertion. 
Now consider equation (23). Again using Stirling’s formula (18) we 


can write 
ioe ae af NS 
M, (y) at eam) (7) : (24) 


If we put y = a - # in (24) (a being a real number greater than 1 and in- 
dependent of /), we obtain: 
ie 
Mr (ah) = ae p" ’ 
where 
p= <1 if «a>t. 


Since for p < 1 the quantity /(h/ 2a) p* tends to zero as h > ©, we must 


have n < ah, or 


i<? for hoo. (23) 
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Inequality (25) can be made to hold for any value of a > 1, no matter 


how near to 1, provided only is chosen large enough. Comparing with 
(20) or (22) we conclude that: 


u] 
h (ho) 1. (26) 


If hk be considered as a continuous variable we can use these results to 
study the derivative dy(h)/dh. Putting: 


OH, (y) _ 


pl 
5) few (1 —n) Indu = Ky (9) >0 (27) 


OH; (y) = fev —w) In (1 = 0) du=L,(v) <0, (28) 
0 


Oh 
we have: 
CS AG Mpuee og 2G), 
"(cd keener mo 
Since — In(1 — u) > u (for 0 < u < 1), it follows that 
yi >1. (30) 
On the other hand, integration by parts yields the formula: 
1 
Hyan (9) = —S 45H) (31) 
or, replacing h by (h — 1), 
ee 
A, (y) = Se ee OPS PETER (32) 
If y = n, so that H,(n) = 1, we obtain from (31) and (32): 
Hy) eee (33) 
n 
and 
tad) 
Ay-1(n) = (=e (34) 


Now, since (for0< u< 1)—(1—u)-M(1—u)<1—-—(1—~»), 
we have: | 


— En (a) < Haas (9) — Hy (9) = ERED (35) 
Furthermore: 
Ky (n) = Hy (0) — Haas (n) =A E DY (36) 
Hence, from (29), (35), and (36): 
fey <itaht2) 7 (37) 


(n—h+1) (h-1) (ee) 


‘IGNITION’ IN A RANDOM NET 373 
This shows that: 


(i>) dh aoe 
It is clear that all the relations found so far will be satisfied if we assume 
that y(/) be of the form: 
y(h) =h+e(h), 
where e(/) is a function tending to infinity with h less rapidly than h/ it- 


self, i.e., 
PE RAED Fie os 


but 


e(h) _ , 
pl, aye ae 0 and also im, Th 


Going back to equation (23) and always assuming h to be a large num- 
ber, we put: 


v=ht+e (39) 
and 
Zao (Ki). (40) 


Then, using Stirling’s approximation once again, (23) becomes: 


ee (1 to] *=V, (41) 
or, putting : 
1 h 
c=Vt in (5-) («1) (42) 
and expanding the left-hand side of (41): 
e~€ [1 tao] b= etl G20? +08) A /)ott—. od, (43) 
o2 o® og C2 
ieee Ae oe 44 
TROPEE Sco ao eye tamey: on 


If in (44) we neglect all powers of o higher than o”, we find: 
oC. (45) 


<) &Vi in (4). (46) 


Going one step further we obtain the equation 


which means: 


Lg eek (47) 


see oc 
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whose solution (correct to the order c?) is: 
ox orton (48) 


Therefore a better approximation for ¢(i) will be: 


= NP (se) +3 In (;-). (49) 


For large h’s the second term in the right-hand side of (49), though still 
becoming infinite with h, is nevertheless much smaller than s/h \n(h/ 2). 
The third term in (49) would be of the order # - c* and therefore negligible. 

Finally, having replaced equation (7) by (23) and comparing the latter 
with (6) we see that A(k) x y(h), or: 


DUDE EEN (se) +3 In (=-) (h>>1). (50) 


Equation (50) gives a much more precise relation than Rapoport’s 
crude estimate 


A (h) > W220 (h—1) 


[Joc. cit., (28)]. Unfortunately, however, it is not a very good approxima- 
tion for values of k which are only moderately large (and which are likely 
to be of greatest interest). Thus for # = 11 the true value of A is 16.93, 
whereas equation (50) gives A = 13.67 (relative error ~20%).* 

Nevertheless, the foregoing considerations are sufficient to show the 
general behavior of the function A(h). As a main result we find: The igni- 
tion phenomenon will not take place unless the axon density is larger than 
the individual threshold (although A and h are of the same order of mag- 
nitude). Also if @ = A and ignition occurs, the fraction y = y(h)/A(h) 
of neurons firing tends to 1 for large values of h. 
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The applicability of the theory of random nets to the theory of rumor spread is 
shown. In particular the ‘“‘weak connectivity’’ of the net appears as the saturation 
fraction of ‘‘knowers’’ in a thoroughly mixed population through which a message 
diffuses, where each knower tells the message to a finite average number of indi- 
viduals. Further it is shown how the time course equation of rumor spread, where 
time is measured by the number of ‘‘removes’’ from the starters, can be translated 
into an ordinary continuous time course equation if the distribution of the telling 
intervals is known. 


R. Solomonoff and A. Rapoport (1951) have defined a random net as 
an aggregate of points connected by directed line segments, such that the 
probability that a connecting link exists between a given pair of points is 
the same for all pairs. They have defined the “weak connectivity” of the 
net as the probability that a chain of connecting links exists between an 
arbitrary pair of points. If the total number of points approaches infinity, 
and if the average number of links issuing from each point is a, then, as 
Solomonoff and Rapoport have shown, the weak connectivity of the net 
is given by the single non-negative root y of the transcendental equation 


2 = tee 8? (1) 


The method used by Solomonoff and Rapoport consisted of “tracing” 
the links issuing from an arbitrary point to the expected number of points 
to which they led, then tracing the links issuing from these to the expected 
number of newly contacted points, etc. i. 
Recent experiments conducted by the Washington Public Opinion 
Laboratory (1952) on the spread of a message through a population fol- 
lowed exactly the same procedure. The points of departure were the 
“starters” of the message. From them the connecting links led to the 
“knowers once removed,” from these again to the “knowers twice re- 
moved,” etc. If the population is “thoroughly mixed,” Le., if the proba- 


SLO 


376 ANATOL RAPOPORT AND LIONEL I. REBHUN 


bility of transmitting the message from one individual to another is the 
same for each pair of individuals, one of whom is a “knower” and the 
other a “non-knower,” then the random net is a suitable mathematical 
model for a theory of rumor spread, provided that the time is measured 
not in ordinary units (hours, days, etc.) but by the number of removes 
from the original knowers. 

In the paper by Solomonoff and Rapoport it is tacitly assumed that 
there is exactly one starter (original knower). As will be presently shown, 
their result is easily generalizable to the cases where the number of start- 
ers is arbitrary. If one speaks of fractions of the population instead of num- 
bers of individuals, then the assumption made by Solomonoff and Rapo- 
port is that the fraction of starters is arbitrarily small. We shall generalize 
to the case in which this fraction is finite. 

Of interest in the theory of rumor spread is the time course of the rumor 
which in the random net model appears as the expected number of points 
contacted in all at each remove (i.e., at each step in the tracing procedure). 
This time course by removes was given by Rapoport (1951), namely, 


Putt =[1- Pi] ue, (2) 


where P(#) is the fraction of newly contacted points. Then the total frac- 
tion of points contacted will be given by 


t 
x(t) = >> P()) (3) 
7=0 
and the expected saturation fraction (ultimately contacted) by 
(0) = S'P(j) =. (4) 
0 


To test equation (2) in any field experiment, it is necessary actually to 
trace the rumor through the knowers at each remove. This was done in the 
C-ville experiment as described by Stuart C. Dodd (1952). In applying 
the data of the experiment to equation (2), two corrections had to be 
made. First, since in the tracing procedure only the actual knowers of the 
rumor were traced, the number of knowers had to be multiplied by the 
average fraction of knowers who became tellers. Second, the quantity a 
(the average number of hearers who became tellers) was not constant 
throughout the process but declined steadily at each remove. 
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With these corrections taken into account and with the actually ob- 
served (corrected) values of a substituted into each recursion indicated by 
equation (2), the comparison shown in Table I between predicted and ob- 
served numbers of tellers at each remove was obtained. 

If, on the other hand, instead of taking the corrected values of a as 


TABLE I 
. Tellers Tellers 
peenieye (Observed) | (Predicted) 

USES oo, ae 42 (starters) 42 
1 ee eee 111 125 
Rarer. Reeve 164 173 
Sac ake ea eee 178 184 
LEASES A ange tiie 180 185 
RR hay 184 186 


calculated from actual observations, one assumes a simple decreasing func- 
tion of a to give the best fit for equation (2), namely, 


Peta 


eens ie (5) 


a 


one obtains a better correspondence between observed and corrected num- 
bers of tellers, as shown in Table IT. 


TABLE II 

Tellers 

Remove ne [By formulas 

(2) and (5)] 
Oe car nce, ees 42 42.0 
i ere ee it ae 111 111.6 
Di wach, Tate e 164 162.8 
Swett nes 178 179.9 
Ee Rade Oe Oe eee 180 183.5 
Swart he at si 184 184.1 


Without further experimentation, nothing can be said concerning either 
the significance of formula (5) or the excellent agreement between theory 
and experiment which it gives. 

Our task in this paper is to derive a method for translating the time 
course by removes formula into an ordinary (continuous) time course for- 
mula. Obviously, the work in tracing a rumor through all the knowers, 
hearers, and tellers is an enormous one, and without it formula (2) is of no 
use in an experimental situation. Before we proceed to derive the transla- 
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tion of (2) into an ordinary clock time course equation, we shall give an- 
other expression somewhat more general than (2) and simpler. Starting 
with equation (11) of Solomonoff and Rapoport Joc. cit., we have 


ay(t) 
y +1) (1-3) ~eamebaak == Ke (6) 


Here (f) is the number of non-knowers at the ‘th remove. We introduce 
a change of notation, in which y will now stand for the fraction of non- 
knowers, and K for the appropriately modified constant. This gives 


aypat n(i-z) = K. (7) 


For very large NV, this can be written as 
y(t+1) eo wH=K. (8) 


To evaluate K we note that if «(é) is the fraction of knowers, so that y(é) = 
1 — x(t), we have 


¥(0)=i—-2; why —= ae (9) 
where «9 is the initial fraction of knowers. Thus 
K= (1—2) e*. (10) 
Equation (7) then becomes, in terms of x(#), 
x (¢+1) =1— (1—-2)) e #4, (11) 


which is simpler and more general than Rapoport’s formula (2). It is an 
immediate consequence of (11) that y must satisfy equation (1) if a_ = 0. 

Equation (11) holds only if ¢ is independent of ¢. If a is a function of é, 
the more general formula (2) still holds. 

We will now show that if the distribution of the time intervals between 
“hearings” and “‘tellings” is known, the clock time course of the rumor 
spread can be expressed in terms of the (7) and the parameters of the dis- 
tribution. 

Call the interval between each time an individual has heard the rumor 
for the first time and the time he has told it (not necessarily for the first 
time) a “telling interval.” Suppose that we have a list of all the telling in- 
tervals in the process studied. Further, suppose that g(r) is a function such 
that g(r)dr represents the probability that a telling interval will have a 
length between 7 and + + dr. Thus g(r) is a probability density and the 
corresponding (accumulative) probability distribution 


G () =f dr (12) 
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represents the probability that a telling has occurred by the time ¢ follow- 
ing a hearing. If r(1) is the fraction of the population who first learned of 
the rumor in the first remove, we see that r(1)G() represents the fraction 
of that fraction completed by the time ¢. We now seek the probability dis- 
tribution of the lengths of chains of telling intervals. If G,(é) is the proba- 
bility distribution of the lengths of chains composed of » telling intervals, 
then we have, under certain assumptions to be stated below, 


Gali =G 7. Gn (13) 


where the operator on ¢ is the m-fold convolution of G, defined by the 
iterative operation 


G, (i) Wea GAC (14) 


Then, if 7(7) is the fraction of the population who become knowers on the 
mth remove, and f(#) is the fraction who are knowers by the time ¢, we 
must have 


f@) = dor (n)G,) = 20+ >) [x (nm) — 2 (n—-1)1G,, (15) 
n=0 n=1 


which is the ordinary time course formula desired. 

The justification of formula (13) depends on the assumptions that the 
individuals act independently of each other and that the tellings of each 
individual are independent. Then the telling intervals will be formed in- 
dependently and according to some distribution G(¢) which will remain 
constant for all removes. The distributions G,(¢) are then the probability 
distributions of the sum of m independent random variables, each having 
the probability distribution G(?). 

Another way of looking at it is to consider the number of knowers JV, 
on the wth remove. Since chains of telling intervals are traced from each 
knower through the first individual who told this knower and so on back 
to the starters, we see that each knower has a unique chain of telling inter- 
vals connecting him to a starter. Thus J, is also the number of such chains 
of m intervals. Hence N,G,(é) is the expected number of such chains of 
length ¢ or less. By the one to one relation between these chains and the 
knowers, V,G,(t) is also the number of knowers on the mth remove at © 
time ¢. If WV is the total population, we have for the fraction of knowers 


at time ¢ 


So MiG, ) Wt = Dr (mG D (16) 
n=0 n=0 


and the total fraction of knowers by the time ¢ is thus given by (15). 
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We will now derive the time course equation for a special case,* namely, 


where : 
g(r) =ke*. (17) 
Then 
; ‘bebdeel (18) 
r(@i) = aa e* i i 
Gi) =a) J 
Similarly, 
Gs (2) = [G (i— 1) dG (1) =Gi@ — Bie, (19) 
0 
G:) =G® —SE (20) 


An easy induction reveals that 

Gr) =6,() —20 eH = 1 — “HE, (Rd), (21) 
where £,,(R#) is the partial expansion of e*' to the mth power inclusive. We 
further define Z_,,(k#) = 0 for 2 > 0 and note that E,)(&#) = 1. Thus also 


Go(t) = 1. Our equation (15) reduces for the special case of the decaying 
exponential frequency distribution of telling intervals to 


f@ = >> r(n)[1— e-*E,_1(kd)1. (22) 


n=0 


But 
rm) =y7, 
n=0 
where y is the saturation frequency as above. Thus we can write (22) as 
f(t) =y— er (n) Ey (he). (23) 
n=0 
We shall now investigate some of the properties of this function. We have 


f' () = ket)” ry [Ena (Rt) — Ena (RE) 
n=1 
: (24) 
= ket fs (Rt) *-1 
n=1 


(n—1)!° 


: Recent data obtained at the Washington Public Opinion Laboratory (unpublished) 
indicate that actual distributions of telling intervals resemble the form we are considering. 
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Clearly f’(¢) = 0 for all? = 0, as, of course, should be the case. Further- 
more, 


ae olan Son we -)). es) 


Consider the last expression in the bracket. This is positive for0 <t< 
1/k. On the other hand, for any mo there exists a fo sufficiently large such 
that 


ki 
(1 -—) =.0 for NSN and BE thy c (26) 


' Since 


= (kt) "71 kt 
Po Oi (: a 


n=1 


is a convergent series, there exists an m such that 


sD Trl at oe (1- =| > P »2, Tati Ae Sa (1- =|. 2.1) 


If we choose 7» so large that 79 > m, then for the value of t) mentioned 
above and for all greater values of ¢, we shall have 


> Trt er (4 =) <0. (28) 


Thus this expression increases at first, reaches a maximum and decreases. 
Let it reach the maximum value M att. Thenif M < 7, we see that f” (#) 
is always negative. If, however, M > nm, then there exist two points ¢’ 
and #’” such that t’ < 4 < ¢’” and f’’(é) < Oforé < #’ and foré > ¢”, and 
f’'() = 0 for t’ < ¢ <¢”. The graph of such a function f(#) is shown in 
Figure 1. 

Note that if 7; is sufficiently small, the initial period of negative second 
derivative is negligible and we have practically a sigmoid curve. On the 
other hand, if 7; is sufficiently large, f (¢) loses its inflection points and be- 
comes a “fall-away”’ curve. In fact, since the expression within the bracket 


in equation (25) can be written as 


ye ee as (r5— rj-1) ; 


we see that we shall always have a “‘fall-away”’ time course curve under 
our present assumption on g(r) ifm > rm Sra... etc., ie., if the ordinal 
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time course curve traced through the data obtained by removes is also a 


“fall-away”’ curve. 
If g(r) is a normal frequency distribution, with mean m and standard 


deviation g, i.e., 


g (7) Se. el— (7— m)?/207} ‘ (2 9) 


o Vln 
the successive convolutions will be given by 


d 1 2 2 
G; (t) ayk. Be el-(r—m)? 20") J 7 


; ; 1 
= ees 5 [—(r¥—nm)?,2no?] dr , 
GW) ie Gay iain s 


f (t) 


f(t)=% 


t 


Ficure 1. Solid line represents the time course of rumor spread for the case 71 > M1, where 


‘ = (ki) "1 ( ) 
Mee Tye, oan ore pr 
with respect to ¢. Dotted line represents the case r1 < MV. 
The integrals (30), being generally available in tabulated form,* the 
values of f(¢) can easily be obtained in this case from equation (15), if the 


r,s are known. 
We note further that 


fO=YV reo. (31) 


*The contributions to the integrals from negative values of ¢ are meaningless but also 
negligible if the mean m is sufficiently large and o not too large. 
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Again f’ > 0 for all ¢, as should be the case. Also 


1 oy tm (nm — t) 
WPAN p= = ¢_  p[—(t—nm)?/2ne?)} 
f) a igo Vn ox 2no2 — ces 


Here, regardless of the values of r, (provided they are convergent), 
f’'(é) = 0 for ¢ sufficiently small and f’’(#) < 0 for ¢ sufficiently large. The 
question of whether and how /” changes sign in the intermediate region 
remains open. 

This investigation is part of the work done under Contract No. AF 
19(122)-161 between the U.S. Air Force Cambridge Research Labora- 
- tories and the University of Chicago. 
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THE DONNAN EQUILIBRIUM 


J. H. Bartrett anp R. A. Kromuout 
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URBANA, ILLINOIS 


A complete solution has been obtained for the distribution of ions and for the 
electrostatic potential in a simple Donnan equilibrium, namely, where the mem- 
brane is permeable to two ions, but not to the third. Electrical neutrality is shown 
to be an extremely good approximation everywhere except close to the membrane. 


If a system contains a membrane which is permeable to some ions but 
not to others, it will finally settle down into an equilibrium state where the 
relative concentrations of the ions on the two sides are markedly different 


membrane 


Ficure 1. System with semi-permeable membrane 


from what they would be if no membrane were present. The ions which 
cannot go through the membrane exert electrical forces on those which 
can, and it is this that causes the difference. 

In order that equilibrium be achieved, let us suppose that the membrane 
(which is imagined to be a geometrical plane, i.e., of zero thickness) spans 
the cross-section of a vessel as shown in Figure 1. This vessel is closed on 
one side, and open to the atmosphere on the other. Since, at equilibrium, 
the concentrations are different, the level of the solution will not be the 


same on one side as on the other. 
385 
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For simplicity, suppose that there are just three types of ion present, 
namely, R+, H+, and Cl-, and that the membrane is permeable to Ht and 
Cl- but not to R+. (One could take into account the difference of osmotic 
pressure on the two sides of the membrane [Guggenheim, 1949], but the 
solution will be regarded as so dilute that this effect is negligible, and also 
so dilute that activity can be replaced by concentration.) Initially, a 
known amount of RClis put in solution on one side of the membrane, and 
a known amount of HCl is put in anywhere. Our problem is then to find 
the equilibrium concentrations of the ions as a function of the distance. 

In this paper we shall show that this problem has a solution which satis- 
fies the following conditions: (a) the electrostatic potential V and its first 
derivative are continuous everywhere; (b) the electric intensity — (dV /dx) 
is zero at the walls x = —L, +L; (c) the electrochemical potential energy, 
per mole, of each ion is a constant; and (d) Poisson’s equation is obeyed. 

The mathematical expression of condition (c) for the ith ion type is 


RT |n ¢;-+ 2,F V = Const. , (1) 


where 
R = gas constant per mole . 


T = absolute temperature , 


c; = concentration of 7th ion , 7 = 1 refers to Ht 
4=2toCl- ,7=3 toR* 


2; = valence of ith ion , 
F = Faraday . 
Equation (1) may also be written as 
Cy = Cyae(-aFV/RT) | (2) 


where c,, is the concentration when V = 0. Let the position of the mem- 
brane be x = 0, and let V = Oat x = —L. 
Poisson’s equation is, if D = dielectric constant, 


a?V 4rF S 
(aa) = 7 Cay De Bee (G)) 


Substitution of (2) into (3) yields 


ad? V 
ia (e2¥ —Je-4V) , 
where 
ae An FP 
oars C20 , Kl = ——(C1a-F C20) 5 coraieg: (4) 
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On integration, we have 


EV Va ae = 
Tin) te (S) 
If we require condition (b) that (dV/dx) = 0 at x = —L, then, from 
(5S), A = —/ — 1. The quantity 7 — 1 may be positive or negative; let us 
study the case / > 1. Then, from (4), at « = —L, (dV /dx”) < 0. This 
means that V < 0 near « = — L. From (4), we see that (d?V/dx?) < 0as 


long as V < 0 and/ > 1. Hence V is always negative and bending away 
from the x-axis, at least up to the membrane. After the membrane is 
crossed, there is no longer any Rt ion, which means that / must be replaced 
by a different constant J). Then, as we shall show, V can bend back toward 
the «-axis, to become asymptotic to a negative value Vo at x = +L. The 
charge density p = —(D/4rm) (d?V /dx?) will then be positive for x < 0 and 
negative forx > 0. The total charge is fp dx = —(D/4r)A(dV/dx). This 
will be zero if dV /dx equals zero at « = —Z and « = +L, and the solu- 
tion will be electrically neutral on the whole. 

Let us now integrate (5) for the general case where (dV /dx) = 0 when 
V = Vo. We set 7? = e*”, and 72 = e*”°. Then we find: 


he = n+(4). (6) 


Equation (5), on rearrangement, becomes 


ate iN eee meas 
(7) 
ae dt 


Now the elliptic integral of the first kind is (Jahnke-Emde, 1938) 


@ dy sing ar 
eRe f aa ary J, Va Ad—eo) 


Integration of (7) gives us 


i 
7 
&— y= + V— =z rou (sin 2, 


and, therefore, 


Z Oly % — Xp xr 8 
¥—Vo= 42 in sn [tv Til: (8) 
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The sn function is 0 at « = x and again at 


yeh 2-0 2K =F (G, z). 
2 2 


T0 


In order to keep V finite, we must have 
Kal) X — Xp ie be oo (9) 
Vig con a2 f VI—P sin’ @ 
Let ho = V1 — 2%. 
From (9), if L = length of the region, 


F 4AcF 
DRT ED. 


As an extreme case, assume L = 10-4 cm, and ¢1, + ¢3, = 10-* mol/ 
liter = 10-7 mol/cc. The constants are taken to be 


D=805 ~F 27.395 Xi eeu, 


2 


DA Sy | saz (C1a+ Ca) i . (10) 


R=8,315 X10? erg! deg|mol, T = 300°. 
Substituting in (10), we find 
Kr eile (11) 


First, consider the region to the left of the membrane, and use equation 
(8), with 


To=1 and k=—= Ti: 


For a value of K as great as 8 [inequality (11)], it is seen from tables that 


kis nearly unity (Jahnke-Emde, Joc. cit., p. 85). In the region (Jahnke- 
Emde, p. 73), 
+ 1 
ee ———) -—-12-— =16. 2) 
Keine; h-layee 

Thus, the value of J) is very close to unity, which means that the solu- 
tion is practically electrically neutral at V = 0. The same argument ap- 
plied to equation (8) for ro = 1 demands that 72~/Jp = 1. The concentra- 
tion of positive ions at Vy is measured by J)/r3 and the concentration of 
negative ions by 74, so that these two concentrations are again very near- 
ly equal at V = V4. If we put Jy) = 1 in (8), then, if r) = 1, 


Vert “in tanh | + NEE (x — 20) | 
(12) 
= Seeintet: ees y, 
a a 


ifw = +V«a/2 (~ — x) and y = tanh xu, 
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The region within which the potential varies turns out (see below) to be 
of the order of 300 A.U. in length. Since Z is assumed greater than 10-4 
cm, it may in this approximation be regarded as infinite. 

Designate the regions to the right and left of the membrane with the 
subscripts 7 and /, respectively. Let us try to fit smoothly two solutions 
of type (12), namely, 
ie 


2 
V ==In tanh | — 5 


(2—a,) |, pxa0h (13) 


V—Vo= —2intann | +V%2(e— 20 ], x>0. (14) 


The first solution vanishes at « = — o, and the second solution van- 
ishes at x = +. Set 


Kya 
uy = vee (41 — x) 


and 
w= VEE (em). 
a a ! du; 
“dx (= q Coth ar sech? ai 7" 
rece Sr «<0, 
ae 15 
zy = 2 cothm, secht u, Fe (15) 


SEN OG ser a eer 
= ee a 


From (13), Vi = (2/a)In 7 and V, — Vo= —(2/a) In yy. 
Equating V; and V, at x = 0, 


Vo== ny, or = Vir = To- (16) 
= = —4es.u. 2 26 mv.) 
(Note: 1/a = RT/F = .862 X 10° e.s.u. = 
Also equating dV/dx|, and dV /dx|, at « = 0 and using (14) and (16), 
i vi ed yi 7 ’ 


1+ 75 _ 1.02136 


vi=— Tae 1200) 


if Vo be taken as —100 m.v. for definiteness. 
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We now find: 


pl Aes 2045: ers ei) 
ie ENT. ies ae 
Up = 201, 


V (eo = 0) = 52a 


Thus, the value of V at x = 0 is —17.47 m.v. if Vo = —100 m.v. 

Concentrations. Suppose that at V = 0, [CI-|-. = 0.1 mols/liter. Then 
at V = —100 mvj [Cla = 0:1) e1" = 2136 < 10 Becauseno: 
electrical neutrality, [H+]+. = 2.136 X 10-%. 


@ 
se 
ae 
sO 
Om 
+ 
= Or 
Do a 
= oO « 
s £, 
= IS0E6 
3 Sy 
= Pa 
o 2 
— 120= 
a o 
Oo 
<7 
°o 
ro) 


Distance From Membrane (Angstroms) 


Ficure 2. Potential and concentrations vs. distance 


At V=0, [H*]-« = 2.136 X 10- XK e108 = 4.56 X 104. Again 
due to electrical neutrality, [Rt]_.. = 0.1 — (4.56 X 10>). 
For 2 
x<0,« =19.5(.897 —u)A, 


x >0, *= 133.2 (w—.2045) A. 


Let 6, be that value of « such that: 


—100—Y¥ (6,) ==[-100—V(u,)]. (17) 
Similarly for 6;: 
0- V6) =210- Vw) 1, (18) 
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or 


V (62) == (-11.17) = —6.31 mv., 


V (6,) = —100+30.45 = —69.55 mv. 


From (17) and (18), the potential thus decreases to 1/e of its value in 
59.4 A if x > 0 and in 9.6 A if <0. 
We may then define the “double-layer” thickness 6 at the membrane: 


6= 6,1 6,=694. 


Thus the “double-layer” contains 28.26 X 6.023 1073 & 1.7 X 1018 
ions/cm’, or a cell of volume 6° in this “double-layer” contains about 10 
ions on the average. This means that any use of a macroscopic theory for 
nonequilibrium calculations—where long time averages of C, V, and their 
derivatives cannot be assumed—would be of doubtful validity. 

Figure 2 shows how the concentrations of the various ions, and the po- 
tential V, vary as a function of distance. This presents then the solution 
of the Donnan equilibrium problem. 
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AWARDS FOR POST-DOCTORAL STUDY IN STATISTICS 
AT THE UNIVERSITY OF CHICAGO 


Three $4000 post-doctoral fellowships in Statistics are offered for 1953- 
54 by the University of Chicago. The purpose of these fellowships, which 
are open to holders of the doctor’s degree or its equivalent in research ac- 
complishment, is to acquaint established research workers in the biologi- 
cal, physical, and social sciences with the crucial role of modern statistical 
analysis in the planning of experiments and other investigative programs, 
and in the analysis of empirical data. The development of the field of Sta- 
tistics has been so rapid that most current research falls far short of at- 
tainable standards, and these fellowships (which represent the third year 
of a five-year program supported by The Rockefeller Foundation) are in- 
tended to help reduce the lag by giving statistical training to scientists 
whose primary interests are in substantive fields rather than in Statistics 
itself. The closing date for applications is February 1, 1953; instructions 
for applying may be obtained from the Committee on Statistics, Univer- 
sity of Chicago, Chicago 37. 
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